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HAUPTAUFSÄTZE 


Über die unter einem Damm durch eine horizontale Parallel- 
schicht sickernde Wassermenge und die Auftriebsdruck- 
verteilung an der Dammbasis. 

Von H. Roßbach in Karlsruhe. 


Das ebene stationäre Problem der Grundwasserbewegung unter einem Damm 
mit Spundwand auf einer horizontalen homogenen Parallelschicht wird 
exakt gelöst. Die Formel für die unter dem Damm sekundlich versickernde 
Wassermenge wird für verschiedene Verhältnisse der Dammbreite, Spund- 
wandlänge und Schichtdicke ausgewertet und in einem Diagramm dar- 
gestellt. 


Einleitende Bemerkungen. In einer Arbeit behandelt T. Kano!') folgendes Grundwasser- 
strömungsproblem: Auf einer homogenen, durchlässigen Parallelschicht von der Breite b ist 
ein Damm aufgesetzt. Die Pegelhöhen des aufgestauten Oberwassers und des Unterwassers 
sind H, und H, (H, — H,—=h). Es wird angenommen, daß die Pegelhöhen konstant gehalten 
werden, und daß der Damm senkrecht zu seinem Querschnitt sehr lang ist; dann kann der 
Strömungsvorgang als ebenes stationäres Problem behandelt werden. Für einen Damm (Basis- 
breite 2a) mit einer und ohne eine senkrechte Spundwand (Länge m) in der Mitte der Damm- 
basis leitet T. Kano Formeln für die Berechnung der unter dem Damm durch die Sicker- 
schicht sekundlich durchsickernden Wassermenge und der Auftriebsdruckverteilung an der 
Dammbasis ab. Ihrer Herleitung nach handelt es sich um Näherungsformeln, die ihre Brauch- 
barkeit verlieren. wenn das Verhältnis der Dammbasisbreite (2a) und der Spundwandlänge (m) 
zur Schichtbreite (b) nicht hinreichend klein ist. Für das vorliegende Problem lassen sich 
jedoch strenge Lösungen angeben, die bei numerischen Rechnungen an Handlichkeit den 
Kanoschen Formeln nicht nachstehen; außerdem können sie auch verwendet werden, um 
den Gültigkeitsbereich der Kanoschen Formeln abzugrenzen. Die Lösung für den Damm 
ohne Spundwand (m=0) (Bild 1b) und außerdem auch die Lösung für eine Spundwand 
allein (a=0) (Bild lc) ergeben sich als Spezialfälle der Lösung für den Damm mit Spund- 


1) Kano, Tokutaro: The Quantity of Water Pereolating under a Dam, and the Distribution of Uplift Pressure 
on it, eonsidering the Depth of the Impervious Stratum., Jap. Journ. Astr. Geoph. Bd. 17 (1939), S. 173 bis 184. 
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b 
Al 
Bild la. Bild 1b. Bild le. 


wand (Bild 1a). Das durch numerische Auswertung der Lösung gewonnene Diagramm (Bild 7) 
gestattet, die durch die Schicht sickernde Wassermenge für beliebige a, b, m und h zu er- 
mitteln. Mit unwesentlichen Abweichungen werden in der vorliegenden Arbeit die von Kano 
verwendeten Bezeichnungen benutzt. 


Bezeichnungen. 
Es bedeuten: 
x,y die Koordinaten eines Punktes des Strömungsbereiches S (Bild la, b, e); 
z=x-+iy eine komplexe Veränderliche; 
u,» die Komponenten des Geschwindigkeitsvektors v der Strömung; 
o©=u-—iv die komplexe Geschwindigkeitsfunktion der Strömung; 
p(x,y) das Geschwindigkeitspotential der Strömung; 
y(z2,9) die Stromfunktion, die längs Stromlinien konstante Werte annimmt; 
die komplexe Potentialfunktion der Strömung; 
p den hydrodynamischen Druck: 
g die in Richtung der negativen y-Achse wirkende Erdbeschleunigung; 
o die Dichte des Wassers; 
v die für das homogen vorausgesetzte Erdreich konstante Widerstandszahl: 


K= den Durchlässigkeitskoeffizienten; 


ovv die pro Volumeneinheit vom Erdreich auf die Flüssigkeit ausgeübte 
Widerstandskraft; 
0 die zwischen zwei Querschnitten im Abstand der Längeneinheit durch 
den Strömungsbereich sickernde Wassermenge. 


Die Differentialgleichungen. Unter der Annahme, daß die Flüssigkeit inkompressibel ist, 
und daß die Geschwindigkeitskomponenten und ihre Ableitungen derart klein sind, da& 
Glieder höherer Ordnung in ihnen vernachlässigt werden können, gelten für den ebenen, 
stationären Strömungsvorgang die Differentialgleichungen 


Aus (la) folgt, daß = als Gradient einer skalaren Ortsfunktion g (x, y) dargestellt werden 

kann, die nach (1b) der Potentialgleichung 
genügt. (1b) kann durch den Ansatz 

befriedigt werden. Wegen 
_dp __dy 


ist 9+iy=f(z) eine analytische Funktion, und es gilt 

= 
Die Lösung des Strömungsproblemes erfordert die Bestimmung einer im Bereich S analytischen 
Funktion f(z), die gewissen durch die hydrodynamischen und geometrischen Daten vor- 
geschriebenen Randbedingungen genügen muß. 

Durch Integration von (1a) erhält man 


y=0 . . . . . . . . . . (2). 


C ist eine absolute Integrationskonstante. (2) dient zur Bestimmung des Druckes. 
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Die Randbedingungen. Längs AB und ÜD ist der Druck » gleich dem statischen Druck 
des anstehenden Wassers. Wegen y=0 ergibt sich aus (2): 


längs AB 


längs CD gen=C- H,. 


Da p nur bis auf eine belanglose, additive Konstante bestimmt ist, kann in einem Punkte 
oder auf einer Aquipotentiallinie willkürlich über 9 verfügt werden. 
Wir setzen daher 


dann wird € 


und R 

Die Dammbasis und Spundwand BEFGC und die Grenze AD der undurchlässigen 
Schicht sind Stromlinien. Auf Stromlinien nimmt y konstante Werte an. 

Die Differenz ypce — yın gibt die zwiscn BEFGC und AD durch den Be- 
reich siekernde Wassermenge Q an. Da auch y nur bis auf eine belanglose, additive Kon- 
stante bestimmt ist, kann in einem Punkte oder auf einer Stromlinie willkürlich über y ver- 
fügt werden. Wir setzen 

und erhalten 

Daraus folgt, daß dem Strömungsbereich S in der @,y-Ebene der Bereich S’ (Bild 2) 
entspricht. Wegen der Symmetrie des Strömungsbereiches zur y-Achse liegt F’ auf der 
y-Achse, E’ und @’ liegen auf der Geraden y'—= Q symmetrisch zur y-Achse. Die Größe @ 
und die Lage der Bildpunkte D’ und F’ sind zunächst unbekannt, sie sind bei der Lösung 
der Aufgabe mitzubestimmen. 


Die Bestimmung der komplexen Potentialfunktion. Die Lösung des Problems verlangt die 
Bestimmung der komplexen Potentialfunktion f(z) aus den angegebenen Randbedingungen; das 
heißt f(z) ist diejenige analytische Funktion, die den Bereich S unter Berücksichtigung der 
Ränderzuordnungsvorschrift auf den Bereich S’ konform abbildet. Im vorliegenden Falle läßt 
sich diese Aufgabe in mehreren Schritten leicht lösen. 


Bild 4. 


Bild 5 (links). 
Bild 6 (rechts). 
G 
1. Schritt: 


2, =e 


bildet den Bereich S auf den Bereich S, (Bild 3) ab. 
Den Punkten A, B, C, D, E, F und @ ordnet sie die Punkte 


gh 
98=5,H,-H)=- 
9 
v | 
G 7 
Bild 2. Bild 3. 
5 
ar m 
zu. 
A 
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2. Schritt: 


bildet S auf S, (Bild 4) ab. Man hat 


itg Tg, 

ig. +? 


= 


| 


3. Schritt: 
(5) 


bildet S auf S, (Bild 5) ab. Die Abbildung liefert derjenige Zweig der Funktion, dessen 
Imaginärteil für alle z aus S größer oder gleich Null ist. Man hat 


4. Schritt: 


(6) 
cos £ 14 Cotg’ tg 

bildet schließlich den Bereich S auf die obere o-Halbebene (Bild 6) ab. In (6) ist für die 
Wurzel das positive Vorzeichen zu wählen. Man hat 


a 
\ 


—1, 


cos Yı + Gig? 


u 


Das elliptische Integral 


do 
\Y 
bildet die obere o-Halbebene auf ein Rechteck der y, y-Ebene ab. Da reellem o((<o< 1) 


0< Ss Eh, y=0 entspricht, muß die Konstante © negativ reell gewählt werden. 


Setzt man (= 2 (e positiv reell), so erhält man 


4293 
itg 395 
2,2, =1, 2, =1, 
| 
a 
Verst; | 
| 
| 
1 2 1 
YııM 
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— k?o®)' 
Da dem Punkte D(o„=1) er y=0 entspricht, erhält man 
1 
do 
— a) (1 
0 


und daraus” h 


(K(k) ist das vollständige elliptische Integral erster Gattung.) Für = erhält man nach 
c 
1 cos’ 


Für die komplexe Potentialfunktion erhält man somit: 


leichter Umformung: 


Cotg,, 


do 
\ (1 — k?o®) 


1 
Da dem Punkte +) q y= entspricht, ergibt sich aus (8) 
1 


do 


1 
1 
\ do +4 \ do 
2» Kik)) Yıl—or)(1— 
0 1 


Nun ist 
1 


k 
\ do 


1 


Damit erhält man für die Wassermenge Q 


An der Dammbasis und den hat man und y=@, daher wird hier 


o 


ya 
k 

1 

und \ =iKık’) 


h. 
69 
I 
yi-P=—. 
a 
ie Cof 
085 
K a 
g Ci 
/ cos’ 
K 1—- 
1 
1) Wegen 
. 
yYııM 
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erhält man 


do 


1 
Im Punkte F(o=x) hat man 9=0. Setzt man 0 = kg 5 erhält man 


1 + Cotg’ 


und 


\ 


Längs FG hat man 


\ 


Längs GC hat man 


3927, +18? 


q 1 ds 
\ 

Aus Symmetriegründen unterscheiden sich die g-Werte in Punkten von BE und EF 
nur durch das Vorzeichen von den p-Werten in den symmetrisch gelegenen Punkten auf FG 
und @C. Unter Benutzung der Druckgleichung (2) erhält man aus den g-Werten die Druck- 
verteilung an der Dammbasis und den Spundwandseiten. Es ist 


h 


Die Spezialfälle (Bild 1b, Bild 1c). Setzt man in den Formeln (8) und (9) u=0 bzw. 
a=0, so erhält man die entsprechenden Formeln für den Damm ohne Spundwand (Bild 1b) 
bzw. für das Wehr (Bild le). Für den Sonderfall «—=0 erhält man aus (8) und (9) 


1 \ do 


1 
a 
- 5, 


2K (20 3) 
Für den Sonderfall «=0 erhält man: 
& 


eos, | 1+ Etg? 3 


do 


und 


mit k=sin 


| 
| 
| 
tg? + 
| 
mit und 
K (eos) 
v u 
(sin 5) 
| XUM 
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Zahlentafel 


Tabelle für die a - Werte. 


0,8 09 


1,0298 0,8071 0,6746 0,5882 0,5000 0,4324 0,3705 0,3097 0,2427 
1.0309 0,9172 0,7724 0,6586 0,5691 0,4952 0,4289 0,3682 0,3082 | 0,2418 
0,8124 0,7763 0,6995 0,6188 0.5453 0.4792 0,4188 0,3613 0,3038 0,2391 
0,6864 | 0,6690 0,6254 0,5701 0,5131 0,4574 0,4037 0,3509 0.2966 0,2343 
0,5991 0,5888 0,611 0,5223 0,4783 0,4324 0,3856 0,3380 0,2877 0,2295 
0,5332 0,5263 0,5071 0,4787 0,4443 0,1064 0,3660 0,3237 0,2775 0,2232 
0,4810 | 0,4761 | 0,4619 0,4402 0,4128 0,3812 0,3466 0,3090 | 0,2671 | 0,2164 
0,4385 0,4347 0,4238 0,4065 0,3842 | 0,3576 | 0,3277 | 0,2943 0,2564 | 0,2094 
0,4124 0.4000 0,3902 0,3773 0,3586 0,3360 0,3099 0,2803 0,2458 0,2023 
0,3728 | 0,3695 0,3631 0,3514 0,3357 0,3163 0,2934 0,2667 0,2356 0,1953 
0.3469 | 0,3449 | 0,3388 0.3288 0,3153 0,2984 | 0,2782 0,2544 | 0,2259 0,1888 
0,2575 | 0.2565 0,2533 | 0,2481 0,2406 0,2309 | 0,2190 0,2042 | 0,1844 | 0,1595 
0.2048 | 0,2041 | 0,2022 0,1988 0,1922 | 0,1798 0,1698 | 0,1567 | 0,1338 


Die Formel (9) wurde numerisch 7 
I, 
für die Werte 
a 1 
und 3 10 
ausgewertet. Die Werte für @ findet 
man in Zahlentafel 1. Nach ihr ist 249 
das Diagramm (Bild 7) hergestellt 
worden. An Hand der Zahlentafel 1 


kann für Werte 0 <1 und 
97 
5=2 und beliebige h, 
und E (Flußbreite) die sekundlich L__2%=05 
versickernde Wassermenge berechnet 
werden. Führt man die Rechnungen 26-07 
für die folgenden Zahlenwerte durch, =08° 
die auch Kano seinen beiden nume- = 
rischen Beispielen zugrunde gelegt hat, | 
Dammbasisbreite: 
2a=60m, Schichtbreite: b= 15m, L__a/b=20 
Pegeldifferenz:: | 
h= 4m, Flußbreitee E=1W0m, 
Durchlässigkeitskoeffizient: | 
K = 0,002 m? sec” ', 
Spundwandlänge: 
m=4m, Bild 7. 


| 


0< 


4 05 


so erhält man für die sekundlich versickernde Wassermenge im Falle m =0: 
0: E=0,02 - 4 100 : 0,2048 = 0,1638 m?. 


Der von Kano errechnete Wert 0: E=0,2244 liegt um 37° über dem nach der exakten 
Formel berechneten Wert. 


Im Falle m =4m erhält man: 
E = 0,2 4 - 100 0,2001 = 0,16008 m?, 


Der von Kano errechnete Wert @- E=0,166 liegt um 3,6°/o über dem nach der exakten’ 
Formel berechneten Wert. 339 


1. 
g h 
v 
m > | 
and 0 0,1 021 08 0.4 0,5 0,6 0,7 1,0 
b 
0 0 \ 
0.1 0 
0,2 0 
0,3 0 
0,4 0 
0,5 
0,6 0 
0,7 v 
0,8 0 
0.9 0 
1.0 0 
1,5 
2.0 0 
| | 
| | | | | | | ] 
| | 
r 
) | | 
N 
— 
\ 
). 
6 070 
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Der Kompressibilitätseinfluß für dünne wenig gekrümmte 
Profile bei Unterschallgeschwindigkeit. 
Von W. Hantzsche und H. Wendt in Braunschweig. 
(Aus dem Institut für Gasdynamik der Luftfahrtforschungsanstalt Hermann Göring.) 


Für die Berechnung ron höheren Näherungen kompressibler Unterschall- 
strömungen um dünne wenig gekrümmte Profile wird ein Verfahren an- 
gegeben, das von der Prandtlschen Näherung für schlanke Profile ausgeht 
und für alle unterkritischen Machschen Zahlen gilt. Dazu wird die Strom- 
funktion nach einem Parameter entwickelt, der die Abweichung der Profil- 
form in Dicke, Wölbung und Anstellwinkel von einem geraden nichtange- 
stellten Streckenprofil charakterisiert. Das erste Glied der Reihe stellt die 
ungestörte Parallelströmung dar. Die beiden ersten (Glieder entsprechen 
der Prandtlschen Näherung. Die neue Verbesserung für etwas dickere 
Profile ergibt sich durch Hinzunahme des dritten Gliedes. Die Methode 
wurde in einigen Beispielen für die maxrimale Übergeschwindigkeit und 
den Auftrieb ausgewertet, die beide stärker mit der Machschen Zahl an- 
wachsen als nach der Prandtlschen Formel. 


Die Berechnung ebener, stationärer Potentialströmungen eines Gases im Unterschall- 
gebiet kann nach Janzen (1913) und Rayleigh (1916) durch ein Verfahren schrittweiser 
Näherung vorgenommen werden. Es entspricht dies einer Entwicklung der Potential- bzw. 
Stromfunktion nach der Machschen Zahl. Die erste Näherung ist also die inkompressible 
Strömung. Die bei dem zweiten und dritten Schritt erhaltenen Approximationen gelten dann 
zwar für beliebige Profile und Anstellwinkel, aber nur für kleine Machsche Zahlen. Strö- 
mungen bei höheren Machschen Zahlen wird man praktisch schwer nach diesem Verfahren 
berechnen können, weil dazu eine größere Zahl von Schritten nötig ist. Die Rechenarbeit 
erhöht sich aber mit jedem Schritte beträchtlich. 

Für dünne, schwach angestellte Profile, bei denen die Zusatzgeschwindigkeiten zur An- 
strömgeschwindigkeit klein sind, haben Prandtl (1927 veröff.) und Glauert (1928) eine 
Näherung angegeben, die sich für alle unterkritischen Geschwindigkeiten anwenden läßt. Die 
Differentialgleichung für das Potential bzw. die Stromfunktion nimmt in dieser Näherung eine 
Gestalt an, die durch affine Transformation einer unabhängigen Variablen auf die der inkom- 
pressiblen Strömung zurückgeführt werden kann. Einer Anregung von A. Busemann') folgend 
haben wir eine Methode der sukzessiven Approximation für die Stromfunktion ausgearbeitet, 
die die Prandtl-Glauertsche Näherung als ersten Schritt hat. Sie entspricht einer Ent- 
wicklung der Stromfunktion nach einem Parameter, der die Abweichung des Profils vom nicht- 
angestellten Streckenprofil charakterisiert. 

Wir stellen zunächst ($ 1) die Differentialgleichungen für die einzelnen Schritte des 
Näherungsverfahrens auf. Für die zweite Näherung geben wir dann ($ 2) die Integration der 
Differentialgleichung an. In $3 wird die zweite Näherung für einige Beispiele ausgerechnet. 
Die Profile sind dabei so ausgewählt, daß sie den Kompressibilitätseinfluß der charakteristischen 
Größen eines Profils (Wölbung, Dicke, Anstellwinkel) erkennen lassen. 


$ 1. Die Differentialgleichungen für die einzelnen Näherungen. 


Die Stromfunktion Y(x,,y) einer stationären rotationsfreien kompressiblen zweidimen- 
sionalen Strömung genügt der Differentialgleichung zweiter Ordnung: 


1 
wobei 
uU Pr; . . . . . . . . . . (2) 
die Komponenten der Geschwindigkeit in Richtung der x- bzw. „y-Achse eines recht- 
winkligen Koordinatensystems sind, und o und a Dichte und Schallgeschwindigkeit be- 
deuten. Der Index © bezeichnet den Zustand im Unendlichen, wo Parallelgeschwindigkeit 
(tz =— 1, angenommen werde. Dazu treten noch die aus dem Adiabatengesetze 
und der Bernoullischen Formel folgenden Beziehungen 
1) A. Busemann: Effeet of eompressibility on thin slightly cambered aerofoils at subsonie speeds. British 


association for the advancement of scienee. Dundee Meeting 1939. Programme 8. 56. Der dort angekündigte Vortrag 
fiel wegen des Krieges aus. 
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Bei Unterschallströmungen ist die Gl. (1) von elliptischem Typ. Für sehr dünne Profile 
kann man im Anschluß an Prandtl ihre Integration auf die für eine inkompressible Strömung 
gültige Gleichung 4 Y=0 zurückführen. Vernachlässigt man nämlich in (1) die Quer- 
komponente v® und setzt «= — [’ der Anströmgeschwindigkeit im Unendlichen gleich, 
so erhält man: 

a; 


Diese Gleichung kann man aber durch eine Transformation der unabhängigen Variablen 


U? 


= 


(6) 
auf die Gestalt Y:+W,,=0 bringen, die für die Stromfunktion der inkompressiblen 
Strömung gilt. Transformation (6) bewirkt eine affine Verzerrung des Stromlinienbildes in 
Richtung der y-Achse. Die neuen Stromlinien gehören zu der inkompressiblen Strömung um 
das ebenfalls mitverzerrte Profil. Die Potentiallinien werden bei diesem Vorgang natürlich 
nicht wieder in Potentiallinien übergeführt. Will man dies erreichen, dann muß man analog 
mit der Differentialgleichung der Potentialfunktion verfahren. Dies ist der Weg, den Prandtl 
ursprünglich vorschlug. Dann gehen aber die Stromlinien nicht wieder in Stromlinien über; 
und insbesondere muß man zu den verzerrten Potentiallinien ein neues Profil (als orthogonale 
Trajektorie) suchen. 

Nur für sehr dünne, vorn und hinten zugespitzte Profile wird man (1) durch (5) ersetzen 
können. Wir wollen im folgenden angeben, wie die durch {5) erhaltene Näherung durch 
sukzessive Hinzunahme der nächsthöheren Glieder verbessert werden kann. Dazu denken 
wir uns Y nach einem Parameter 4 entwickelt, der die Abweichung der Profilform von einem 
geraden nichtangestellten Streckenprofil charakterisiert (Dicke, Wölbung, Anstellwinkel): 


Die Anströmgeschwindigkeit ist dabei in Richtung der negativen »-Achse angenommen. 

Aus (7) erhalten wir für die Komponenten der Stromdichte 7 =735 = 

Entwicklungen nach Um die Geschwindigkeitskomponenten und selbst 

zu erhalten, verwenden wir die aus (3) folgende Beziehung zwischen der Dichte 


und —— 
Uo 


und der 


) 
Stromdichte 
2 | D 2 D 2 2 


1 
und P=" +(2 +1) gesetzt ist. Mit (7) ergibt sich dann: 


au 


2 


> 


e 
e 
d 
e 
n 
n 
it 
1- 
ie | 
1- 
d 
t, 
t- | 
3 
t. 
wobei 

dx 1 

— 

) [2 

ish 
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Geht man mit diesen Ausdrücken in die Differentialgleichung (1) ein, so erhält man 
auch hier eine Entwicklung nach 4: 


+ gun 


ar 
3 5 
Ta +(&— 1) a’ 
Da die Differentialgleichung für beliebige 4 richtig sein muß, kann man die Ausdrücke, 


die bei A, 4? usw. als Faktoren stehen, einzeln gleich Null setzen. Man erhält für die ersten 
drei Koeffizienten: 


2a? +1] u? 9, — hr) 


Die Differentialgleichungen (9), (10) und (11) gestatten g, h und f nacheinander zu ermitteln. 
Die Randbedingungen für g, h und f erhalten wir dabei aus den Forderungen über 
Y= const auf dem Profil, 
U im unendlichfernen Punkt, 
c) glatter Abfluß an der Hinterkante des Profils bzw. Staupunkt. 


$ 2. Die erste und zweite Näherung. 


Wir beschränken uns im folgenden auf die Berechnung der ersten und der zweiten 
Näherung, behandeln also die Gl. (9) und (10). Führt man an Stelle von g und h die Ab- 
schnitte 7” und 7 der Stromfunktion Y ein: 


so genügen 7‘ und 7® offenbar den folgenden Gln. (12) und (13) 


(1497) 


Da die abgebrochenen Entwicklungen Y und Y® im allgemeinen die Randbedingungen 
von Y nur noch in ausreichender Näherung erfüllen, empfiehlt sich für die ersten beiden 
Schritte folgender anschaulicherer Weg: 

Wir bestimmen eine im allgemeinen von 7" y) verschiedene „Stromfunktion* (x, 
die der Gl. (12) genügt: 


und dabei die Randbedingungen für Y exakt erfüllt, also: 


a) P® = const auf dem Profil, 


b) =— im Unendlichen, 


c) glatter Abfluß an der Hinterkante bzw. Staupunkt. 


wobei 
_ 
XUM 
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Diese Funktion (x, die offenbar der Prandtlschen Näherung gleichkommt, setzen 
wir auf der rechten Seite von (13) für 7“ ein und bestimmen nun eine Funktion Y® (x, y), 
die jetzt (13) erfüllt: 


xx 


. (13a) 
— a —1) U 


und dabei wieder den Randbedingungen von Y genügt: 
a) Y®=const auf dem Profil, 
b) =0, =— U im Unendlichen, 
ce) glatter Abfluß an der Hinterkante bzw. Staupunkt. 
Es ist zu erwarten, daß die so bestimmten Funktionen Y® und 7 sich von PF® und ?® nur 
in Gliedern höherer Ordnung als / bzw. #? unterscheiden. 7“ genügt ja der gleichen Differential- 
gleichung wie 7°, nur wurden dabei die Randbedingungen genau erfüllt. In Gl. (13a) haben 
wir rechts gegenüber (13) zwar 7 durch 7 ersetzt, eine Änderung der Differentialgleichung 
ist dabei aber erst in Gliedern der Ordnung 4° vorgenommen worden. Der genaue Nachweis, 
daß in den Entwicklungen nach 4: 


ya 


die Funktionen (x, y), 9° (@,y) und g‘” y) identisch sind und weiter (x, y) mit h (x, y) 
zusammenfällt, ist zwar grundsätzlich einfach, soll aber hier, da er etwas umfangreich ist, 
unterdrückt werden. 
Um Gl. (12a) und (13a) zu integrieren, gehen wir durch die Transformation 
Y 


(14) 


in die &, 7-Ebene über. Wir erhalten dann: 


(16). 
— a! u" (x 


In der &,»-Ebene haben wir ein Profil, das durch ähnliche Verzerrung längs der -Achse - 
aus dem Profil der x, y-Ebene entsteht. Die Gl. (15) definiert y” (£,,) als Stromfunktion 
der inkompressiblen Strömung um dieses Profil. Die Anströmgeschwindigkeit dieser inkom- 
pressiblen Strömung ist mit Rücksicht auf und == Y wieder — U. 

Um 7°” (&,n) und y” (£,7)) zu bestimmen, bilden wir weiter das Profil der £&, 7-Ebene 
auf den Kreis r=yr?+9y’=R der komplexen 3-Ebene durch eine konforme Abbildung des 
Außengebietes des Profiles auf das Außengebiet des Kreises r=R ab, wobei der unendlich 
ferne Punkt der &,»-Ebene in den der 3-Ebene übergeht. Die Abbildungsfunktion sei durch 
die komplexe Beziehung gegeben: i 


Mit 


wird nun aus (15) und (16) mit det je 


12 
e, 
), 
), 
). 
n. 2 

2), 
3). 
en 
en 
| 
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(18). 


Für K® bzw. K“ haben wir offenbar die Randbedingungen: 
K®’=const bzw. const auf dem Kreis. 
—() =—TDi J ic 
dh U im Unendlichen. 
=0, im Bildpunkt der Hinter- 
kante des Profiles. 


K“ (x,9) bestimmt sich nun wegen (17) als m U WG), wenn wir mit UWG) das 
komplexe Potential für die inkompressible Strömung um den Kreis der 3-Ebene mit der An- 
strömgeschwindigkeit — 7 bei dem durch die Hinterkante des Profiles bestimmten Staupunkt 
bezeichnen. Mit den Symbolen m bzw. Re wollen wir hier und im folgenden den Imaginär- 
bzw. Realteil einer komplexen Größe bezeichnen. 

Die Integration von (18) läßt sich mit Hilfe von W(3) und der Abbildungsfunktion £ (3) 
ebenfalls angeben. Es gilt zunächst offenbar: 


Wegen de "arde kann man nun für (18) wie man sich leicht überzeugt schreiben: 
- 


/ 


wenn man zur Abkürzung unter $ bzw. y: 


dWd3z 


dl 


versteht. Der Querstrich bedeutet Übergang zur konjugiert komplexen Größe. 

Zur Lösung von (19) empfiehlt es sich an Stelle von x, y die komplexen Variablen 3=r+1Y 
und 5=r—iy einzuführen. Bei der Behandlung der inkompressiblen Strömung hat schon 
Busemann?) diese Transformation verwandt. Imai und Aihari’) haben sie neuerdings 
bei der Durchführung der Rayleighschen Näherung benutzt. 

Mit dieser Transformation nimmt nun (19) die Gestalt an: 


21) 


dt 


ergibt die Integration dann: 
® 2) A. Süsamann in A. Föppl: Vorlesung über Technische Mechanik Bd. IV, Dynamik. 8. Auflage, Leipzig 
1933, S. 383, 
3) ]. Imai und T. Aihari: On the bee flow of a eompressible fluid past an «elliptic eylinder. Rep. of 
Aeron. Research Inst. Tokyo, Bd. XV, 8 (1940), S. 185. 


YııM 
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+ 
Die Funktionen F (3) bzw. G 65) sind in 5 bzw, 5 analytische Funktionen. F(6)+G(6) be- 
stimmt sich durch die Randbedingungen für K” (S. 76). 


Die formale Verwendung von 3 und 5 als unabhängige Variable bedarf noch der Recht- 
fertigung. 


Der konjugiert komplexe Wert f(5) einer analytischen Funktion &) von 5 ist keine 
analytische Funktion von 3 mehr; dagegen ist f(3) offenbar analytisch in 5. Wir schreiben 
um dies auszudrücken ER 


Wegen d;—d;j gilt dann offenbar für die Ableitungen von f nach ;: 
(25) 
Damit kann man Gl]. (19) mit Rücksicht auf (22) schreiben: 
(26). 
Läßt man jetzt De bo ER Variablen x und 9 jede für sich kom- 
plexe Werte 1, +ir, und 9, +9, annehmen, so werden offenbar (5), W6), 6). WG) und 


K“ analytische Funktionen der beiden komplexen Variablen x und y. An Stelle von x und y 
führen wir nun als unabhängige Veränderliche 


und 


ein. Dann geht (26) über in: 


| -( a du d;, -{ d;. "9, da, | 


Für K“ ergibt sich also: 
2 2 4 


wobei F(3,) und @(3,) analytische Funktionen von 3, bzw. 3, sind. Beschränkt man sich 
wieder auf reelles x und y, dann wird 3,=3 und ,=3;. Nimmt man weiter, um eine reelle 
Lösung von (19) zu erhalten, noch den Realteil von (27), so bekommt man schließlich mit 
Rücksicht auf (24) und (25) den Ausdruck (23): 


Wir wollen dafür 


”(2) 


Bei der Bestimmung von F(5)+ @() empfiehlt es sich, zunächst die Funktion W (65) 
| Im w) abzuspalten; dadurch ist die Bedingung im Unendlichen erfüllt. Den Rest 


| 

), 
) 
zig 
of 
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ermitteln wir in der Weise, daß wir eine Summe von Funktionen von 3 bzw. 3, die im Außen- 
gebiet des Kreises r—R bis auf den unendlichfernen Punkt analytisch sind, mit folgenden 
Eigenschaften suchen: Mit entgegengesetztem Vorzeichen soll 

1. der Imaginärteil des Restes mit dem von H (3,3) auf dem Kreise r—= R überein- 
stimmen, 

2. die Ableitung des Imaginärteiles nach x und y im Unendlichen gleich der von 
Sm H (3,5) sein. (Wegen der Bedeutung von ß, W und £ verschwinden übrigens die ersten 
Ableitungen von ‘m H (3,3) im Unendlichen.) 

Durch Hinzunahme eines zu ‘m (ö- In 3) proportionalen Gliedes kann noch der Staupunkt 
in Ordnung gebracht werden. 

In den meisten Beispielen ersetzen wir die Bedingung 1. durch die weitergehende 
Forderung, daß sowohl der Real- als auch der Imaginärteil von H (3,3) auf dem Kreise voll- 
ständig durch die Zusatzfunktion, die man dann leicht angeben kann, kompensiert wird. 
H (5,35) zerfällt nämlich in eine Summe von Produkten analytischer Funktionen der Form 
h,(3)h,(3). Für ein einzelnes Glied h, (3) - h, (3) kann man dann durch Hinzufügen des Ausdruckes 


—h .(@ :) h,(5) bzw. —h, 65) () stets erreichen, daß auf dem Kreise die Randbedingung er- 


füllt ist. Dabei treten aber unter Umständen Singularitäten im Außengebiet des Kreises 
auf. Immer ist dies nicht der Fall, wie durchgeführte Beispiele zeigen, auch ist oft eine 
einfache Beseitigung dieser möglich. 


$ 3. Die Berechnung der zweiten Näherung. 


Nach dem eben entwickelten Verfahren ist der zweite Schritt für den Kompressibilitäts- 
einfluß in einigen Beispielen gerechnet worden. Die Profile sind so ausgewählt, daß sie den 
Einfluß der charakteristischen Größen Wölbung, Dieke und Anstellwinkel zeigen. 


a) Einfluß der Wölbung (Nichtangestelltes Kreisbogenprofil). 
Wir untersuchen zunächst ein nichtangestelltes Kreibogenprofil in einer Parallelströmung. 
Die Abbildungsfunktion 


(28) 


bildet den Kreis r=a der 

j,Ebene_auf einen Kreis- 

bogen AB der Pfeilhöhe f 

in der £ Ebene ab. Die 

Punkte A und B sind die 

Bild Ib. Hilfsebene. Bilder der dureh y=— £ 

festgelegten Punkte A’ und 

B’ des Kreises in der 5 

Ebene. In der x y-Ebene 

erhalten wir nach (14) 

dann einen „Kreisbogen“ 

y 7, B,mitder Pfeilhöhe fu. 

(Bild 1a, b, 

Bild 1a. Kreisebene. Die komplexe Strom- 

5, funktion U W der inkom- 

pressiblen Strömung mit 

der Anströmgeschwindig- 
keit — U um den Kreis r—a in der xy-Ebene bei den Staupunkten A’ und B? ist: 


Bild le. Strömungsebene. 


UW=- Ufilng 


Der Parameter A, die Wölbung des Kreisbogens A, A, B, ist 7= = 4 


Daraus ergibt sich f oder bis zu den quadratischen Gliedern in 


YıımM 
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f= — . Wir haben also in (28) und (29) den Wert f = zu wählen, um in der x, y- 
Ebene um ein Kreisbogenprofil der Wölbung r eine Strömung mit der Anströmgeschwindig- 


keit — U zu erhalten. Der Einfachheit halber entwickeln wir nach f oder £ und schreiben 


erst das Ergebnis auf 7 um. Der Bogen A,B, in der «y-Ebene ist genau genommen ein 
Ellipsenbogen, die Abweichung von einem Kreisbogen der Pfeilhöhe fu liegt aber erst in 
der Größenordnung (f «)’, kann also bei der Berechnung des zweiten Schrittes vernach- 
lässigt werden. 


r 


Die Größe wird nach (29): 


U 
1) 1 a? 
Dabei wurden bei den Termen W)ß bzw. Jm(£ 5°) nur die Entwicklung nach f bis 


zu den Gliedern in f? angegeben. Die Beschränkung darauf ist gestattet, weil unser zweiter 
Schritt nur die Glieder bis einschließlich f? richtig wiedergibt. 


Die Funktionen F(53) und @ (53) lassen sich nach dem S. 78 angegebenen Verfahren so- 
fort angeben. K“” nimmt dann die Gestalt an: 
K® a? 


4 4 +D 


Inz ist dabei dadurch eindeutig gemacht, daß die Ebene längs der negativen x-Achse 
aufgeschnitten ist. 


Damit erhält man — — 


Daß die Punkte 3= + % a? — R 3 Staupunkte sind, bestätigt man leicht. u rg ist 


in diesen nach (29) Null. Die mit f? behafteten Glieder K® geben aber zu e 5 für 


— 
=% Ve- _ f 5 nur Beiträge in höheren Gliedern als f?, denn die unter den Zeichen 
‚m stehenden Ausdrücke sind nach r abgeleitet für 3= + a reell. 


4ra 
Setzt man 3=reif, so erhält man daraus in reeller Form, wenn noch f=- 


geführt wird: 
2 4 
=(- r+ — Inr 
2 72 n? 4 


/ r 


2 2 4 


(2 +1) a) 


16r 16r° 
Aus dem Ausdruck für K“ ergibt sich dann wegen 


und K® —y® bis 


zu den Gliedern 7? genau die Näherung 7” der Stromfunktion. Berechnet man nun daraus 
die Maximalgeschwindigkeit am Kreisbogen, so erhält man für diese einen bis zum quadra- 
tischen Glied in r richtigen Ausdruck: 


| 
| 
5 in; 
| 
| 
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Bild 2 zeigt die Maximalgeschwindigkeit an einem 
Kreisbogen der Wölbung 1 in Abhängigkeit der Mach- 


M Machsche Zahl der \ | schen Zahl. Kurve 1 entspricht dem Ausdruck (30), 
ungestörten Strömung | u] wenn man nur das konstante und lineare Glied in r 
| | berücksichtigt. Kurve 2 und 2a stellen den voll- 
ständigen Ausdruck (30) für =1 bzw. «—=1,4 dar. 
Die daneben noch eingezeichneten Kurven 1’,2’ usw. 
entsprechen der Entwicklung der eigentlichen Prandtl- 
schen Näherung nach r (12)). 
Die Kurven € geben an, wo die kritische Mach- 
sche Zahl für =1 bzw. #=1,4 erreicht wird. 
Der Auftrieb ist in A“ durch die Vorzahl von 
In r bestimmt (vgl. Anhang). Da diese sich gegenüber 
der in K‘” nicht geändert hat, ist der Auftrieb der 
gleiche wie bei der ersten Näherung, in Überein- 
daß er eine ungerade Funktion der 
2 neue Methode (7.bzw 2 Schritt) grıst 
Der Kreisbogen stellt das einfachste Beispiel für 
unsere Methode dar. In diesem -Falle kommt man 
auch leicht mit reellen unabhängigen Variablen x,» 


angestellten Kreisbogenprofil der Wölbung durch. 4 
Gl. (19) nimmt nämlich die Gestalt an: 


ar a "sing 


Bild 2. Maximalgeschwindigkeit bei dem nicht- 


(31). 
+1) „[sin3g  a’sing 


Nun überzeugt man sich aber leicht, daß AK=Ar"’sinng (A,p,n Konstante) das 
folgende partikuläre Integral hat: 
ap*?r 
(P+2)°- (pP —n r" 
K=-Asinng 2(p+2) | 

Da diese Ausdrücke für K auf dem Kreis r—= a konstant sind und im Außengebiet des 
Kreises keine Singularitäten haben, kann man damit ohne weiteres die Lösung von (31) an- 
geben; sie stimmt mit der nach der komplexen Methode erhaltenen überein. 

Bei den folgenden Beispielen ist die komplexe Methode vorzuziehen. Durch Entwick- 
lung der rechten Seite von (19) nach Termen der Form r? cosng bzw. r’ sinng könnte man 
auch mit reellen Variablen arbeiten. Nur ist die Konvergenz der Entwicklung in den Bei- 
spielen nicht immer gewährleistet oder es sind auch nur Entwicklungen nach Doppelsummen 
möglich, so daß man praktisch auf diese Weise schwer durchkommt. 


b) Einfluß der Dicke (Ellipse, Kreiszweieck). 


1. Nichtangestellte Ellipse. 

Wir berechnen nun den zweiten Schritt für die nichtangestellte Ellipse in Parallel- 
strömung. 


K=Asinng wenn (p+ 2? — 


wenn (p 2)? :— 


bildet den Kreis der eine Ellipse mit Achsen: 


a+d— und a+d+— 


der Ö-Ebene ab. Die nachfolgende Verzerrung (14) gibt in der x y-Ebene eine Ellipse mit 
den Achsen: 


a? 
u +d— und a+ 


oder dem Achsenverhältnis 


(a+d)? a? 


07 
DI 065 
950 
055 
050 
045 
035 
er 
) 
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; d 
Wir werden im folgenden d bzw. = als Parameter verwenden und erst das Ergebnis 


auf e umschreiben. 
Die komplexe Stromfunktion U W() der Strömung um den Kreis mit dem Radius 
r—a+d, den Staupunkten 9=0, und der Anströmgeschwindigkeit — ist: 


+ Ad)? 


Mit (32) und (34) wird 
bzw. K“® bekommen nun die Gestalt: 
K' a* (2 + I) a? (a A)? a d rs 
a' u +1). a? ad+.d? 
Die Funktion F (5) + @ (6) läßt sich nach dem Verfahren von S. 78 angeben. Man erhält dann: 


K“ a? u? (a+d} 


= im WG) = — 


(36). 

Man erkennt leicht, daß vr die Stromfunktion für die zirkulationsfreie Strömung um 
die nichtangestellte Ellipse ist. Bei Spiegelung an der x-Achse (Übergang zur konjungiert 
komplexen Größe in den hinter dem Zeichen ‘m stehenden Ausdrücken) wechselt K” nur 
sein Zeichen. 

Im Gegensatz zur Behandlung des Kreisbogens haben wir in (35) keine Entwicklung 
nach dem Parameter d vorgenommen. Das wäre auch hier möglich und würde die Rechnung 
vereinfachen. So nimmt (35) z. B. dann die Gestalt 


( +D ta? 
[7 x > 


an. Dies ist aber nicht empfehlenswert, weil man nun im Rahmen dieser Näherung die 
Randbedingung für die mit d? behafteten Glieder am Kreise r—=a und nicht bei r=a+d 
erfüllen wird. Auf diesem Kreise liegen aber die singulären Punkte der rechten Seite von 

Aus dem Ausdruck (36) für K“® kann man nun die Maximalgeschwindigkeit an der 

Ellipse ermitteln. Mit (33) = = 
hältnisses e bis zum quadratischen Glied ine genau: 


- ergibt sich für eine Ellipse des Dickenver- 


Bild 3 zeigt die Übergeschwindigkeit bei einer Ellipse der Dicke z in Abhängigkeit 


von der Machschen Zahl der Anströmung. Die Kurven © geben wieder an, wo die kritische 
Geschwindigkeit erreicht wird. 

Die Bilder 4a und 4b zeigen noch die kritische Machsche Zahl als Funktion des 
Dickenverhältnisses der Ellipse. Zum Vergleich sind dabei wieder die von Kaplan nach 
Rayleigh-Poggi gerechneten Werte eingetragen. Die weiter eingezeichnete Kurve von 
Tsien*) ist mit der Hodographenmethode gewonnen, wobei die adiabatische Druck-Volumen- 
beziehung linear angenommen wurde (Tangente in dem dem Zustand der Anströmung ent- 
sprechenden Punkte). 


D H. s. Tsien: Two-dimensional subsonie flow of eompressible fluids. Journ. Aeron. Sei. Bd. 6 (1939), S. 399 bis 407. 
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M Mochsche Zahl der ungestörten | 


net ine: 


2Schritl) 


| 


a 8 8 0 


0 


Bild 3. Maximalgeschwindigkeit an einer Ellipse 


mit dem Achsenverhältnis 


T 
beschwindigkeit 


a” kritische Deschwindigkeit 


| 


| 
| 
| 


(peue Methode 1 Schritt) 
| 
AS. Tsien 
| | 
76 
"Dicke 
/iefe 


Die kritische Anströmgeschwindigkeit 2 


a* 


Bild 4a. 


als Funktion des Achsenverhältnisses der Ellipse. 


U Anströmgeschwindigkeit 
a* kritische beschwindigkeit 


| 
Methode 2. Schritt 


N 


1Schrit 
Kaplan (Rayleigh 2.Schritt) 


A357) 


pressiblen 


Bild 7. 


% Bild 4b (links). 
Übergeschwindigkeit 


03 


Aurve 1: neue Methode 1.chritt (Prandil 
„ 2: neue Methode 2 Schritt (n = 14) | 
„ 3: Methode 2.Schritt = 7) | 
Prondil Glied in und | 

| | 


403 


Bild5. Maximalgeschwindigkeit beim niehtangestellten 


Kreiszweieck der Dieke 


Bild 6. 


als Funktion der Dicke des Kreiszweiecks. 


Die kritische Anströmgeschwindigkeit 


Joukowskiprofil, Kreisebene bzw. Strömungsebene. 


Differenz der aus der kompressiblen und inkom- 
berechneten kritischen Anström- 
geschwindigkeit bei der Ellipse. 
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2. Nichtangestelltes Kreiszweieck. 


Den Einfluß der Dicke haben wir noch an dem Beispiel des nichtangestellten Kreis- 
zweieckes ermittelt. Die Abbildungsfunktion, die den Kreis r—=a in der ;-Ebene auf das zur 
&- und y-Achse symmetrische Kreiszweieck in der £-Ebene mit dem Kantenwinkel und der 


(n=2-2) ib durch 


Stna \ta 
gegeben. Das senkrecht zur Anströmrichtung affıin verzerrte Kreiszweieck wurde wieder durch 
ein Kreiszweieck approximiert. Wir können dies mit der gleichen Berechtigung wie beim 
nichtangestellten Kreisbogen tun. Die Rechnung führen wir hier nicht näher aus, sondern 

geben nur das Ergebnis an. Bezeichnen wir mit 7 das Verhältnis Tiefe , 50 erhält man für 
die Maximalgeschwindigkeit am Kreiszweieck bis zu den quadratischen Gliedern in 7 genau: 


4 12 2 5 


In Bild 5 ist für = die Maximalgeschwindigkeit in erster und zweiter Näherung 


angegeben. Die Übergeschwindigkeiten für eine Ellipse mit dem gleichen Diekenverhältnis 
liegen etwas tiefer (Bild 3). Dies ist aber inkompressibel schon der Fall. Vergleicht man 
das Kreiszweieck mit einer Ellipse, deren Dieke so gewählt ist, daß die inkompressiblen Über- 
geschwindigkeiten bis zu den Gliedern in 7’ gleich werden, so zeigt sich auch für höhere 
Machsche Zahlen eine sehr gute Übereinstimmung, ein Ergebnis, das von vornherein nur 
für die Prandtlsche Näherung zu erwarten ist. 

In Bild 6 ist weiter der Eintritt der Schallgeschwindigkeit in Abhängigkeit von der 
Dieke des Kreiszweiecks angegeben. Zum Vergleich ist noch die entsprechende Kurve für 
die dünne nichtangestellte Ellipse eingezeichnet. 


ec) Einfluß des Anstellwinkels (angestelltes Jouko wskiprofil). 


Zur Ermittlung des Anstellwinkeleinflusses haben wir für das angestellte symmetrische 
Joukowskiprofil in Parallelströmung die durch die Kompressibilität bedingte Auftriebs- 
änderung untersucht. Das durch die affine Verzerrung (14) aus einem Joukowskiprofil 
hervorgehende Profil wurde dabei wieder durch ein Joukowskiprofil approximiert. 


Ein angestelltes Joukowskiprofil ist charakterisiert durch die aus Bild 7 ersichtlichen 
Parameter: a (Radius des umströmten Kreises), ö (Anstellwinkel), e (Verschiebung des Mittel- 
punktes des Einheitskreises gegen den des Stromkreises). Um die zweite Näherung für die 
Strömung um dieses Profil zu erhalten, nähern wir das nach (14) affın verzerrte Joukowski- 
profil wieder durch ein solches an, dessen Parameter «’, ö’, e’ seien. Setzt man: 


so stimmt das Näherungsprofil in den linearen Gliedern in e und ö mit dem verzerrten 
Profil überein. (Von der Abweichung in den quadratischen Gliedern könnte u. U. noch das 
Glied in eö in (42) einen Zusatzterm in e geben. Nach (43) ist dieser sicher von mindestens 
der Größenordnung «*.) Wir führen die Berechnung zunächst in a’, ö’,e’ durch und schreiben 
das Ergebnis dann auf a, öd,e um. Da wir hier bei festem a’ zwei Parameter « und 
haben, denken wir uns diese zur Durchführung unseres Verfahrens linear an 4 gekoppelt 
—=&,4, —6,4 und die Stromfunktion nach entwickelt. 


Das komplexe Potential UW für die Strömung um ein Joukowskiprofil der Para- 
meter a’, ö’,e’ bei der Anströmgeschwindigkeit — U ist 
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wobei die Abbildungsfunktion {(3) die Gestalt: 


hat. Aus (38) und (39) ergibt sich 5 zu: 


(je-i® +’) — a’ 
jeri®d a’ +28 
Damit nehmen K‘” und K“ die Werte an: 


K“’ a’? 


“(2 +N 


U S 


(3 


(a eidr 
x 


— a’ +2: 


— a’ +28 


3 eier + ıe 2 


Wählt man zur Erfüllung der Randbedingung AK“ = const auf dem Stromkreise nach 
96) 
- , ?e-id 


dem Verfahren von S. 78 die Funktion von }, so tritt ein Glied 96) ist im Punkte 
m. 5 regulär) und damit ein Pol im Außengebiet des Kreises auf. Durch Hinzufügen von 


wird dieser aber beseitigt. Da wir uns auf den Auftrieb beschränken, genügt uns nach dem 
Anhang die Kenntnis der Vorzahl des Gliedes in Inr der Funktion K“ in der ;-Ebene. Der 
Auftrieb und damit diese Vorzahl ist aber sicher eine ungerade Funktion in ö’, enthält also 
nur Glieder in ö’,e’ö’,.. usw. Wir brauchen also in (41) nur die Glieder in ö° und e’ö' zu 
bestimmen. Die Vorzahl .{ von \\m[itn 5] ergibt sich dann aus der Forderung, daß der 
Punkt 3= — «’ ei” Staupunkt werden soll zu 


a N 


> 


A=—- 2a — + | 


Mit (37) wird daraus: 


au 


2eöl, 1 
7 u u 


Für das Verhältnis der Zirkulation 7/7; der kompressiblen zu /'; der inkompressiblen 
Strömung erhält man dann nach Formel (8) des Anhanges: 


4 1 


also (bis zum linearen Glied in = genau 


- 


= ist bis auf einen Faktor angenähert das Dickenverhältnis: 


Dicke 


| 
| 
| | 
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> 
ui Das Verhältnis r, hat Kaplan’), der nach 


Rayleigh-Poggi (Entwicklung nach der Machschen 
Zahl) vorging, für dünne Profile zu: 


angegeben. Dies ist das erste Glied der Entwick- 
lung (42) nach a. In Bild 5 ist der Ausdruck (42) 
für verschiedene Diekenverhältnisse d in Abhängig- 
keit der Machschen Zahl aufgetragen. Die Kurve 
d=0 entsprieht der Prandtlschen Auftriebsformel: 

N; 


= u. 


Anstell- | 
winkel 


Die Berechnung des zweiten .Schrittes unseres 
Näherungsverfahrens zeigt in allen Beispielen, daß 
der Kompressibilitätseinfluß durch Dicke, Wölbung 
Bild 8. Auftrieb des symmetrisehen und Anstellwinkel gegenüber der Prandtlschen 


Joukowskiprofiles bei kleinem Anstellwinkel 
und kleiner Dicke. Näherung noch etwas ansteigt. 


Anhang. 


Berechnung des Auftriebes aus der Stromfunktion. 


Die Stromfunktion hatte in der zweiten Näherung bei den vorhergehenden Rechnungen 
mit 3=reir für große r eine Gestalt, die sich der Form: 


r 


n, 


unterordnen läßt. Dabei seien «a, und 5b,„, Funktionen von g, außerdem gehen in sie wie 
auch in die Konstante 1 die das Profil charakterisierenden Parameter a, d,.,d,f usw. und 
die Machsche Zahl a ein. m,n und k sind ganze Zahlen = 1. 


Die Funktion wird, da die Abbildung 


(e e,**+* Konstante) das Unendliche unverändert läßt, wieder die Form (1) haben. Mit 


sich: 
k 


n, = 
Die Abbildung = gibt zwischen den Polarkoordinaten r,# der &n-Ebene und 
R, ® der x y-Ebene die Beziehungen: 
sin? 


tg 


Damit nimmt y)= u y(£,n) die Gestalt an: 


m 


Ü 
n, 


wobei die @,„ und b,; Funktionen von ® sind. 
Aus dem Ausdruck (3) berechnen wir nun die Komponenten der Stromdichte. 


man die Glieder durch Punkte an, die im Unendlichen mit höherer Ordnung als RB 
werden, so erhält man: 


6) C. Kaplan: Compressible flow about symmetrieal Joukowski profiles. NACA Rep. 621 (1938), 8. 8. 
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on 


„Asin® a’ u .Icos’ Psin 
Udy i 


+ | 


R 


R cos? ® 


/ 2 \2 
Die Dichte ergibt sich aus der Stromdichte (5 +[ e nach der aus 


Gl. (3) von $ 1 folgenden Entwicklung nach 1 | = ): 


o 
Rlcos®’ PD + — 


ın D 2 D sı 
a’ u Acos’ d sin 1 
R sin?’ 


Nach dem Adiabatengesetz erhält man schließlich für den Druck u =(2 ): 


On 


(6). 


2,72 


p Lcos’ Psin ® 
Ricos’ P-+ 


Legt man in der xy-Ebene einen großen Kreis R=C um das Profil und wendet auf 
die Flüssigkeitsmenge im Kreise den Impulssatz an, so findet man für den Auftrieb P des 
Profiles: 


wobei das Integral über den Kreis R=([ zu erstrecken ist. cos (nz) und cos (ny) sind die 
Richtungskosinus der inneren Normalen, ds das Linienelement auf dem Kreis. 


Die rechte Seite von (7) nimmt nun mit (4), (5) und (6) für R>x die Form an: 


a? cos’ sin’ 


27 
— U’uA P+ 
Mit z=«a’o, a, =o, Ü” wird daraus: 


a 


d 


u” 


was ausgewertet: 
P=-o 


gibt. Setzt man P=—o, UT'(l' Zirkulation), so hat man: 
Me 


Man sieht daraus, daß bei einer Stromfunktion der Form (1) (bzw. (3)) nur die Kon- 
stante bei dem Glied mit Inr in den Auftrieb eingeht. 325 
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Eine Anwendung des absoluten Parallelismus 
auf die Schalentheorie. 


Von F. Reutter in Karlsruhe. 


Nach einer kurzen Entwicklung der Theorie des absoluten Parallelismus 
werden mittels dieses mathematischen Apparates die Gleichgewichtsbedin- 
gungen für das deformierte Schalenelement aufgestellt und anschließend 
insbesondere die Fundamentalgrößen zweiter Ordnung der deformierten 
Schalenmittelfläche errechnet. 


Der absolute Parallelismus ist eine verhältnismäßig junge mathematische Entdeckung 
und hat bisher in erster Linie in der (mehrdimensionalen) höheren Differentialgeometrie 
Anwendung gefunden. Die vorliegende Arbeit hat nun das Ziel, die Grundgleichungen der 
Schalentheorie mittels des mathematischen Apparates des absoluten Parallelisınus herzuleiten. 
Es handelt sich um die Beantwortung folgender Fragen: 

1. Wann ist eine durch ihre Mittelfläche gegebene Schale unter der Einwirkung be- 
stimmter äußerer Kräfte und Momente im Gleichgewicht, d. h. welche Spannungs- 
resultanten und -momente stellen sich ein? 

2. Welche Deformation erleidet die Mittelfläche der Schale? 


Damit Gleichgewicht besteht, muß die Resultierende der an den Seiten eines Flächen- 
elements der Mittelfläche angreifenden Kräfte und Momente einschließlich der äußeren Kräfte 
und Momente der Nullvektor sein. Zwecks Ermittlung der Resultierenden sind diese Kräfte 
und Momente durch eine Parallelverschiebung in denselben Angriffspunkt zu verlegen. Für 
eine Platte ist diese Parallelverschiebung trivial und für eine Schale mit gekrümmter Mittel- 
fläche läßt sie sich im euklidischen Raum durchführen. Dennoch gelangt man, wie wir hier 
zeigen wollen, erst dann zu einer organischen Entwicklung der Grundlagen der Schalentheorie, 
wenn man sich des absoluten Parallelismus (auf der Mittelfläche der Schale) bedient. Auf 
diesem Wege ergeben sich auch die Fundamentalgrößen der deformierten Mittelfläche wesent- 
lich leichter als bei der von Love durchgeführten Herleitung‘) mittels beweglicher Koordi- 
natensysteme auf der Mittelfläche und in größerer Allgemeinheit als bei der von Trefftz?) 


gegebenen Herleitung mittels des Castiglianoschen Prinzips. Der hier angewandte mathe- 
matische Apparat ist zwar komplizierter, läßt aber die wesentlichen Gedankengänge klar 
hervortreten. 

Die Abhandlung zerfällt in drei Teile. Die Theorie des absoluten Parallelismus wird 
im Abschnitt I entwickelt und im Abschnitt II zur Aufstellung der Gleichgewichtsbedingungen 
benutzt. Der Abschnitt III behandelt die Deformation der Mittelfläche. 


I. Der absolute Parallelismus. 


1. Durch einen Punkt P der Ebene sei eine Gerade y gelegt. In der euklidischen 
Geometrie gibt es dann durch einen beliebigen Punkt P’ der Ebene eine und nur eine 
Parallele y’ zu g. Dieser Begriff des Parallelismus kann sofort auf alle Flächen vom 
Krümmungsmaß Null übertragen werden, also — 
außer der Ebene — die Kegel, Zylinder und die 
Tangentenflächen beliebiger Raumkurven, die Torsen. 
Die ©' Erzeugenden einer Torse (Bild 1) sind be- 
kanntlich zugleich die ©' Schnittgeraden ihrer &' 
aufeinanderfolgenden Tangentialebenen; indem man 
nacheinander diese Tangentialebenen um ihre 
Schnittgeraden dreht, wird die Torse auf die Ebene 
abgewickelt. Daher heißt ein Vektor t in der Tangen- 
tialebene des Punktes P einer Torse zueinem Vektor t’ 
in der Tangentialebene des Punktes P’ parallel, Bild ı. 
wenn nach Abwicklung der Fläche in eine Ebene die 
Bilder t und f’ der beiden Vektoren euklidisch parallel sind und wenn | =|t!=|t|=|t|. 

Jetzt sei eine Fläche ® vom Krümmungsmaß K=F0 gegeben. P und P’ seien zwei 
Flächenpunkte. Der Vektor t in der Tangentialebene des Punktes P soll nach dem Funkte P’ 
parallel verschoben werden. Dazu verbindet man P mit P’ durch eine beliebige reguläre oder 
wenigstens stückweise glatte Flächenkurve C. In allen Punkten von C werden die Tangential- 


1) Love: Lehrbuch der Elastizität, deutsch v. A. Timpe, Leipzig und Berlin 1907, $ 326. 
2) Trefftz: Ableitung der Schalenbiegungsgleichungen mit dem Castiglianoschen Prinzip, Z. angew. Math. 
Mech. Bd. 15 (1935), S. 101 bis 108. 
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ebenen an die Fläche ® konstruiert. Die Einhüllende dieser einparametrigen Ebenenschar 
ist eine Torse 7, die die Fläche ® längs € berührt. Der absolute Parallelismus auf 
der Fläche ® wird jetzt nach Levi-Civita®) mit Hilfe des Parallelismus auf der Torse T 
so definiert: 

Ein Vektor t in der Tangentialebene des Flächenpunktes P und ein Vektor t’ in der 
Tangentialebene des Punktes P’ sind auf der Fläche ® in bezug auf dieKurve parallel, 
wenn die beiden Vektoren auf der der Fläche ® längs C umschriebenen Torse T parallel sind. 
Die Lage des Vektors t’ hängt also von der Verschiebungskurve € ab. Wird die Parallel- 
verschiebung einmal längs einer Kurve €,, das andere Mal längs einer Kurve (', vorgenommen, 
so erhält man in P’ i. a. zwei voneinander verschiedene Vektoren t/, t,: sie haben zwar 
gleiche Länge, schließen aber einen von den Krümmungsverhältnissen von ® abhängigen 
Winkel 7 ==0 ein. Denn die der Fläche ® längs der Flächenkurve €, umschriebene Torse T, 
ist mit der ® längs (€, umschriebenen Torse 7, dann und nur dann identisch, wenn ® das 
Krümmungsmaß Null hat; dann gilt sogar T, = T, =®. Auf Flächen vom Krümmungsmaß 
Null ist die Parallelverschiebung vom Verschiebungswege unabhängig. 

2. Wir wollen nun aus den Komponenten des Vektors t bei gegebener Fläche ® und 
Verschiebungskurve C die Komponenten des parallel verschobenen Vektors t’ berechnen. 

Die Fläche ® sei in vektorieller Darstellung 

gegeben, die Flächenkurve € habe die Parameterdarstellung 


Piu,(s),u,(s)) und P’(u, (s+ds),u,(s+ds)) seien zwei 
benachbarte Punkte von und ihre Tangential- 
ebenen mit der Schnittgeraden g (Bild 2). D ist ein 
Vektor in der Richtung von g. Der Vektor t mit dem 
Anfangspunkte P in r soll durch die Parallelverschiebung 
in den Vektor t’ mit dem Anfangspunkte P’ in r’ über- 
gehen. XL, und tu, sind die Tangentenvektoren an die 
Parameterlinien der Fläche im Punkte Pi(w,,u,). Nach 
Bild 2. Nr. 1 soll |t|=|t’| sein oder auch 
Dreht man die Tangentialebene ı’ um g, bis sie mit der Ebene r NEE NON so müssen 
nach Nr. I die Vektoren t und t’ euklidisch parallel sein. Daher gilt weiter: 
Umformung von (1) und (2) liefert: 
®, 
Wegen 
t=1(s)= t, + t, 
+ d,(s) 


.. 
folgt weiter nach Ausführung des Grenzübergangs ds>0 mit t— re 


Dies ist ein System von zwei linearen EEE erster VE für die Kom- 
ponenten t,(s), t,(s) des Vektors t. Unter Benutzung des aus (3) folgenden Wertes für t 
nehmen (5) und (6) die Form an: 


(7) 


d,9,+4,9,=0 


“=l,:=2; 


» Levi-Civita: Nozione di parallelismo, Reudiconti di Palermo, Bd. 42 (1917). (1. Veröffentlichung) — Der 
absolute Differentialkalkül, Deutsch v. A. Duschek, Berlin 1928, S. 35 ff. — Bezüglich der hier gegebenen Ableitung ver- 
gleiche außerdem: Komerell: Theorie der Raumkurven und krummen Flächen, Berlin und Leipzig 1931, Bd. Il, 
u en Zur geodätischen Krümmung und Parallelverschiebung, Jahresber. der Deutschen Math. Vereinig. 

7 (1928), S. 257. 
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Fällt t nicht in die Richtung von d '), so ist 


t, 
d,,d, 


Daher gestattet das System (7) nur die triviale Lösung 
0,=0 
Die Auflösung von (8) nach (s), t,(s) liefert: 
Hier bedeuten 7”, die Christoffelschen Dreiindizessymbole zweiter Art, aus der Diffe- 
rentialgeometrie bekannte Linearkombinationen der metrischen Fundamentalgrößen der 


Fläche ®, E= (tu, Yu), F= (Xu, Kur), (Ku, und ihrer Ableitungen. 
Die Gln. (9) können auclı age werden: 


du; 


Bei gegebener Fläche ® und Verschiebungskurve € sind die (7 2 u) bekannte Größen: die 


(9). 


GlIn. (10) sind daher ein System linearer homogener Differentialgleichungen 1. Ordnung für 
die unbekannten Funktionen f;. (s), aus dem sich zu irgendwelchen Anfangswerten f}. in einem 
Flächenpunkte P(u;(s)) die Komponenten f,(s) des parallelverschobenen Vektors in allen 
Punkten der Verschiebungskurve € bestimmen lassen. Die Gln. (10) gelten für Mannigfaltig- 
keiten M,„ beliebiger Dimensionszahl n. Da die Koeffizienten der Gln. (10) nur von den 
metrischen Fundamentalgrößen der M,„ und der Gleichung der Verschiebungskurve abhängen, 
ist der absolute Parallelismus eine Riehtungsübertragung, die allein durch die Metrik der M,, 
bestimmt ist und in keiner Weise von ihrer Lage in einem Einbettungsraum abhängt. 


II. Die Gleichgewichtsbedingungen. 


3. Eine Schale wird durch zwei gekrümmte Flächen F\,, F, begrenzt. Deren Abstand 2h, 
die Schalendicke, ist gering gegen ihre sonstigen Abmessungen, kann im übrigen aber be- 
liebig veränderlich sein. Diejenige zwischen F,, F, liegende Fläche ®, die überall die Schalen- 
dieke halbiert, wird Mittelfläche genannt. Durch Vorgabe der (undeformierten) Mittelfläche ® 
mit der Vektorgleichung 

r=r(u,v) 
und der Schalendicke 2h für jeden Punkt von ® ist die Schale geometrisch bestimunt. 

Unter der Einwirkung äußerer Kräfte mögen die Punkte der Schale gewisse Verschie- 
en aus ihrer ursprünglichen Lage erfahren. Diese Verschiebungen — “# in der Richtung 
‚ w in Richtung r= v 
dicke 2h so klein sein, daß die Produkte der Verschiebungsgrößen , ?, w und ihrer partiellen 
örtlichen Ableitungen vernachlässigt werden können. Zur Kennzeichnung des Spannungs- 
zustandes der Schale dienen statt der auf die Flächeneinheit eines Schnittes bezogenen 
Spannungen die Resultanten der über die ganze Schalendicke verteilten Spannungen und 
Spannungsmomente. 

Es sei € eine Kurve auf der Mitteltläche der deformierten Schale mit der Parameter- 
darstellung u = u (s), (s), t ihr Tangentenvektor im Punkte Pis), » die Flächennormale 
im Punkte P mit |t'= »/=1 und n=[t»v]. Durch t werde ein Normalschnitt der deformierten 
Mittelfläche gelegt und in der Normalschnittebene ein Flächen- 
stück begrenzt durch » und den Normalenvektor », im Punkte P, 
der Normalschnittkurve o (Bild 3). Die auf das schraffierte 
Flächenstück in Bild 3 wirkende Spannung läßt sich bei hin- 
reichend kleinem Kurvenbogen PP,=do durch eine Kraft und 
ein Kräftepaar in P(s) darstellen; beide werden in Komponenten 
nach dem System t, ı, » zerlegt. Division durch do und der 
Grenzübergang P,>P, also do —>0 liefert die auf die Kurve € 
bezogenen Spannungsresultanten und Spannungsmomente, die 
auch mit Hilfe der Komponenten des Spannungstensors durch 
Integration über die Dicke der Schale dargestellt werden können 


4) Fällt t in die Richtung von d, so ist 1" mit t schon euklidisch parallel. In der Grenze geht die Riehtung von d 
in die zur Tangentenriehtung an die Verse hiebungskurve „konjugierte Riehtung“ über. 


" in der Richtung sollen gegen die Schalen- 
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(Nr. 4, Gln. 11); der Grenzwert für die »-Komponente des Kräftepaares verschwindet immer. 
Statt der Bezeichnungen: „Spannungsresultante“* und „Spannungsmoment* werden in der 
Literatur vielfach die Bezeichnungen: „Schnittkraft“ und „Schnittmoment“ gebraucht. 


4. Die Parameterlinien der undeformierten Mittelfläche seien Krümmungslinien. Diese 
bilden ein orthogonales Netz, so daß F=V (Nr. 2). Die undeformierte Mittelfläche ist daher 
durch Angabe von E, @ und der Fundamentalgrößen zweiter Ordnung 


Ku u Yu X.) (iv Lu 


L= 
charakterisiert. 


Wir nehmen an, daß die Dehnung der Mittelfläche bei der elastischen Deformation der 
Schale so klein ist, daß sie bei der Ermittlung der Gleichgewiehtsbedingungen vernachlässigt 
werden kann. Dann geht das Netz der Krümmungslinien infolge der Deformation in ein 
Kurvennetz über, das mit hinreichender Annäherung weiterhin als ein Rechtwinkelnetz an- 
gesehen werden kann, aber wegen der bei der Deformation eintretenden Verbiegung nicht 
mehr aus Krümmungslinien besteht. Die deformierte Mittelfläche mit der Vektordarstellung 


kann daher angenähert durch die alten metrischen Fundamentalgrößen E, G, aber neue 
Fundamentalgrößen zweiter Ordnung 


E= (Zuu Yu X.) M= Eu v .) 


charakterisiert werden. Die Spannungsresultanten und Spannungsmomente werden jetzt auf 
das System 

Lu - [Eu 

yE’ y6' ’TYEG 


bezogen. Im Punkte P(u,v) der Mittelfläche seien die zur v»-Linie gehörenden Spannungs- 
resultanten mit 
die zur «-Linie gehörenden mit 
—8,T,, N, 


bezeichnet, die Spannungsmomente mit 


H,, 6, bzw. —G,, 
Ist 


+h 


5) ]. a. ist noch 03 =0 (vgl. Nr. 13, Gl. 21). 


| 
Tıa Tız 
Ta Ta 0) 
der Spannungstensor, bezogen auf das System JR ve. »°), so gelten die Gleichungen : 
—ıh 
2 (r 
—h 
h E + h h E u). 
‚ z-R z:E 
—h 
+ h E +h 
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5. Die Punkte P(w,») und R(u+du,v+dr) 
der deformierten Mittelfläche bestimmen durch ihre 
«- und v-Linien ein kleines Parallelogramm, gebildet 
aus den Kurvenstücken Cu+au, von der 
Länge Y@dv bzw. yEdu (Bild 4). Es besteht Gleich- 
gewicht für das durch dieses Parallelogramm darge- 
stellte Schalenelement, wenn die Resultierende der auf 
diese Kurvenstücke wirkenden Kräfte bzw. Momente #furdu,) u+du Rfurdu,v+dv) 
und der äußeren Kräfte bzw. Momente der Nullvektor ist. Bild 4. 


Es wirken die Kräfte 
mit dem Angriffspunkt P(u,r) auf das Kurvenstück C,. 
8, = (s: | Ic +N, du 
mit dem Angriffspunkt P (u,v) auf das Kurvenstück 


= (r, / | u, d 


v@e (N, y@a)d v(u+du, dr 
mit dem Angriffspunkt P,(u+du,v) auf das Kurvenstück 


N) 
+ vE+ (N,yE)d 0) r(u,o+dr)du 


mit dem Angriffspunkt P, (u,v+d») auf das Kurvenstück ı»- 


Die Orientierung ist so gewählt, daß 8, 8,, S,, 0, aus dem Schalenelement PP,RP, 
herausweisen. Die zugehörigen Momente sind: 


= E-\ 
(6:5 


Zur Ermittlung der Resultierenden müssen die Vektoren 


bzw. M, vom Punkte P,(u+du,v) längs C, 
nach dem Punkte P(w, v), 

die Vektoren $, bzw. M, vom Punkte P,(u,e+dr) längs C,„ 
nach dem Punkte Pu, v) 


gemäß Abschnitt I parallel verschoben werden; ihre Komponenten sind damit auf das 
LyKvv-System des Punktes P(u,v) bezogen. Die äußeren Kräfte und ihre Momente 
lassen sich bei hinreichend kleinem d«, d» durch eine Kraft und ein Kräftepaar mit dem 
Angriffspunkt in P(u,v) darstellen. Damit Gleichgewicht am Schalenelement PP,RP, be- 
steht, müssen die Komponenten der beiden Resultierenden einzeln verschwinden. Durch 
Grenzübergang du >0, dv>0 erhält man schließlich die Gleichgewichtsbedingungen für den 
Punkt v). 

6. Durch die Einführung der Spannungsresultanten und -momente sind zwar alle in der 
Schale wirkenden Kräfte und Momente durch Vektoren mit Anfangspunkt auf der Mittel- 
fläche dargestellt, doch sind dies keine flächenmäßigen, d. h. in der Tangentialebene der 
Mittelfläche liegenden Vektoren. Daher muß zur Ausführung der in Nr. 5 geforderten Parallel- 
verschiebung eine dreidimensionale Mannigfaltigkeit M, konstruiert werden. Ist X (u,») der 
Ortsvektor der deformierten Mittelfläche, so hat diese M, die vektorielle Darstellung: 
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+ w =r(u,r) 


Hieraus erhält man nach bekannten Formeln der Ditferentialgeometrie: 
r 
und die metrischen Fundamentalgrößen 4;, dieser M, mit 


Die Verschiebungskurven €,, ©, (Nr. 4) liegen auf der Mittelfläche, die aus den 9;, zu be- 
stimmenden /-Symbole (Nr. 2) werden also nur in den Punkten der Mittelfläche benötigt und 
in den Differentialgleichungen der Parallelverschiebung (Nr. 2, GIn. (10)) ist daher de =du, 0. 
Hat der zu verschiebende Vektor t im System » die Komponenten t,, so lauten 
die Differentialgleichungen der Parallelverschiebung: 


1 
—op(Eu du+E,de)t, du —G,„de)t,+ (Ldu+Mdv)t,, 


(Ldu+Mdv)t, — (Mdu+Ndo)t, 


Die Vektorkomponenten #; sind der Reihe nach mit den Komponenten von X, $,, 
M,, % zu identifizieren und für die Parallelverschiebung längs €, ist noch dv —0, längs (€, 
ist du=0. Daher ergibt sich z. B. für die Verschiebung von $, längs C..: 


E,, 


/ G 


Die Glieder | au usf., die beim Grenzübergang du >0 bzw. wegfallen, 


sind in diesen schon fortgelassen. 
Nach Ausführung der Parallelverschiebung gehen 
die Vektoren $,, M, in P, in die Vektoren $,, M, in P, 
die Vektoren $,, M, in P, in die Vektoren $,, M, in P_über. 
Ist das Parallelogramm PP,RP, hinreichend klein, so lassen sich die auf seine Fläche 


wirkenden äußeren Kräfte ersetzen durch eine Kraft X in P(w,v) und ein Kräftepaar Q. 
Die Kraft X hat die Komponenten: 


X,y@Gdude, X,yEdudr, X,yEGdudr bezüglich I», 


Damit Gleichgewicht am Schalenelemente PP,RP, besteht, müssen die Komponenten der 
Vektorsumme 
HR, +X 
einzeln verschwinden. 
Bei der Aufstellung der Gleichgewichtsbedingungen für die Momente ist zu beachten, 
daß die auf zwei gegenüberliegende Seiten unseres Elementarparallelogramms wirkenden 
Kräfte in ihren entgegengesetzt gleichen Anteilen Anlaß zu den zusätzlichen Momenten 


geben. Ist Q das Moment der auf die Fläche PP,RP, wirkenden äußeren Kräfte mit den 
Komponenten 


y@dude, Q,yEdudv, 0 bezüglich r, 


cn 
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so müssen daher auch die Komponenten von 
HM 

einzeln Null werden. 

Die Komponente von K in der Richtung r,. wird 


2E E 2E 
NyE+X v6 dude. 
Daher ist A, die 1. Gleichgewichtsbedingung für das Parallelogramm PP,RP,. Durch 
Grenzübergang du >0, dv>0V, — Glieder erster und höherer Ordnung mit du, dr in der 
[]-Klammer sind bei den Rechnungen in Nr. 5, 6, 7 von vornherein weggelassen — erhält 
man schließlich nach leichter Umformung ingesamt die N 6 Gleichgewichtsbedingungen: 
du E 


VE), Er +O)YEG=0, 


du dv ayE 


8. Die 6 Gln. (13), (14)°), zu denen noch für jedes vorliegende Einzelproblem die Rand- 
bedingungen’) angegeben werden müssen, sind 5 lineare partielle Differentialgleichungen erster 
Ordnung und eine gewöhnliche lineare Gleichung für die 10 unbekannten Funktionen 7,,T,, 
N,, @,,@,, H,, H,; das Problem ist also statisch unbestimmt. Bei vorgegebener 
Mittelfläche sind zwar die Größen E, 6, L, N bekannt; die Größen L, M, N lassen sich aber 
nach Abschnitt III nur bei bekannten Verschiebungen der Punkte der Mittelfläche angeben. 
(Da insbesondere an „dünnen“ Schalen auch bei kleinen Deformationen größere Krümmungs- 
änderungen auftreten können, würde die Aufstellung der Gleichgewichtsbedingungen in den 
Größen L, N, also für das nichtverbogene Schalenelement, u. U. zu erheblichen Ungenauig- 
keiten Anlaß geben.) Wie bei einem allgemeinen elastischen Problem aus den Gleichgewichts- 
bedingungen 

X = x-Komponente der äußeren Kräfte 
allein die Spannungen nicht errechnet werden können, so sind mittels der Gln. (13), (14) 
allein die Spannungsresultanten und -momente nicht zu ermitteln. Hierzu muß bekannt- 
lich das Elastizitätsgesetz der Schale herangezogen werden. Mittels der Gln. (11) sind die 
Spannungen einzuführen, diese aber vermöge des Hookeschen Gesetzes und der besonderen 
Annahmen über den Deformationszustand der Schale — Abschnitt III — schließlich auf die 
Verschiebungen der Punkte der Mittelfläche zurückzuführen‘). Ebenso sind die Größen 
L, M, N nach Abschnitt III durch die Verschiebungen auszudrücken, wobei die in Nr. 3 ge- 
machten Annahmen über die Größenordnung der Verschiebungen nach dem Einsetzen in die 
Gln. (13), (14) noch zu Vereinfachungen Anlaß geben. So erhält man endlich ein System 
partieller Differentialgleichungen für die Verschiebungen v), (u, v), (u, v), das freilich 
i. a. überhaupt nicht oder doch nur näherungsweise integriert werden kann. — Bei den 
meisten bisher in der Literatur behandelten Problemen der Schalentheorie ergeben sich aber 
durch die spezielle Gestalt der Mittelfläche, geeignet veränderliche Dieke der Schale usf. so 


6) Die G!n. (13), (14) stimmen mit den bei L or, Ioe. e., S. 612/13 auf anderem Wege hergeleiteten Gleichungen 
überein, wenn man die dort auftretenden Größen ihrer üächenmäßigen Bedeutung entsprechend durch 
E,G,L, M, N ausdrückt. 

Love, 1. e. $ 332. 

s) In der unter ?) zitierten Arbeit von Treiftz, S. 108, sind die Spannungsresultanten und -momente dureh die 
Verschiebungen ausgedrückt, allerdings ausgehend von den am undeformierten Schalenelement aufgestellten Gleich- 
gewichtsbedingungen. 
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weitgehende Vereinfachungen im Ansatz, daß der hier in voller Allgemeinheit angegebene 
Weg nicht oder nur teilweise beschritten werden muß. So reduziert sich bekanntlich bei 
biegungsfreier Beanspruchung einer Schale die Anzahl der Spannungsresultanten — unter 
Wegfall der Momente — auf drei und man erhält ebensoviele Gleichgewichtsbedingungen. 
Derartige Aufgaben sind also statisch bestimmt (Membranspannungszustand)). 

9. Entsprechend dem in der linearen Elastizitätstheorie üblichen Verfahren und den 
Annahmen über die Größenordnung der Verschiebungen (Nr. 3) wurde bei der Aufstellung 
der Gleichgewichtsbedingungen gefordert, daß die durch die Deformation bedingte Dehnung 
der Seiten des Rechtecks PP,RP, vernachlässigt werden kann. Insbesondere soll der rechte 
Winkel zwischen den Koordinatenlinien sieh nur so wenig ändern, daß durchweg mit hin- 
reichender Annäherung die metrischen Fundamentalgrößen E, @ benutzt werden dürfen. Trifft 
diese vereinfachende Annahme nicht mehr zu, so treten in den Gleichgewichtsbedingungen 
Zusatzglieder auf, die auch die metrische Fundamentalgröße F’=(f„X,) enthalten. Diese 
Zusatzglieder lassen sich wieder gemäß Nr. 6 und 7 berechnen. In den Gln. (13), (14) sind 
ferner überall E, @ durch die entsprechenden Größen E’, G’ der verzerrten Mittelfläche zu 
ersetzen. Das in der ersten Gleichung des Systems (13) auftretende Zusatzglied hat z. B. 
die Form: 


+3E’'@ EN, +3MyE@@' N, +NE” N.| 


Jetzt sind nicht nur die Spannungsresultanten und -momente und die Größen Z, M, N, sondern 
auch die vorläufig ebenfalls unbekannten Größen E’, F’, @’ gemäß Nr. 10, Gl. (18) durch die 
bekannten Größen E, G, L, N und die Verschiebungskomponenten x, v, w auszudrücken, um 
ein System von Differentialgleichungen in den letzteren zu erhalten. Außerdem ist zu be- 


achten, daß die Spannungsresultanten und -momente jetzt auf ein stem 7 PR . 


ve 
rechter mehr. 


Lu Iv 


ist aber kein 


‚» bezogen sind. Der Winkel zwischen den Richtungen 


II. Die Deformation der Mittelfläche. 


10. Die Mittelfläche der undeformierten Schale wurde durch die Größen E, G, L, N 
charakterisiert, die Mittelfläche der deformierten Schale durch die Größen E’, F’, @, L,M, 


Die Verhältnisse =, s kennzeichnen die infolge der Deformation eintretende Dehnung der 


Linienelemente der Mittelfläche, F’ die Winkeländerung der Parameterlinien, die Größen L, 
M, N geben gegenüber Z, N die Krümmungsänderungen der Parameterlinien und damit die 
Krümmungsänderung der Mittelfläche an. Der Deformationszustand der Schale ist daher be- 
stimmt, wenn die Größen E’, F’, @', L, M, N durch E, @, L, N und die Verschiebungskompo- 
nenten «, v, w ausgedrückt sind. 

Zur Herleitung der Gleichgewichtsbedingungen (13) und (14) waren entsprechend der 
Tatsache, daß aus ihnen allein die Spannungsresultanten und -momente i. a. nicht ermittelt 
werden können, besondere Annahmen über die Art der Deformation der Schale nicht erforderlich. 
Die GIn. (13), (14) gelten daher auch für Probleme der Schalentheorie, bei denen die nach- 
stehenden Annahmen nicht erfüllt sind. 

Für die Deformation der Schale soll das folgende physikalische Modell gelten: 

Die Linienelemente, die zur undeformierten Mittelfläche normal sind, sollen gerade 
bleiben und zur deformierten Mittelfläche normal werden. Diese Voraussetzung ist immer 
erfüllt, wenn die Schalendicke 2% klein ist gegenüber den Abmessungen der Mittelfläche "). 
Hierzu tritt gewöhnlich die — für unsere weiteren Entwicklungen nicht wesentliche — An- 
nahme, daß diese Linienelemente keine Dehnung erfahren sollen. 

Ist ve) der Ortsvektor der undeformierten Mittelfläche, 

der Verschiebungsvektor, so lassen sich die Punkte der deformierten Schale — bzw. der in 
Nr. 6 eingeführten M, — auch durch den Ortsvektor 


%) Flügge: Statik u. Dynamik der Schalen, Berlin 1934, Abschn. II bis IV; Lagally: 7. angew. Math. Mech. 
Bd. 4 (1924), S. 377; Löbell: Z. angew. Math. Mech. Bd. 7 (197), S. 463. 
10, Flügge: |. e. S. 112. 
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beschreiben. Die Parameterlinien und die Normalen der undeformierten Mittelläche be- 
stimmen durch die Vektorgleichung 

ein krummliniges orthogonales räumliches «-r-w-Koordinatensystem. Die Flächen w = const 
sind Parallelflächen zur Mittelfläche, die Flächen « = const und ®e=const bilden zwei Scharen 
abwickelbarer Regelflächen: deren Erzeugende sind die Normalen der undeformierten BR 
fläche in den Punkten der verschiedenen Krümmungslinien. Sind allgemein a;(i=1,2,3) 
orthogonale Koordinaten eines räumlichen Bereichs, g;, die zugehörigen metrischen Funda- 
mentalgrößen und 9; die zur Tangentenrichtung der «;-Linie gehörigen Verschiebungs- 
komponenten, so erhält man nach bekannten Formeln '') für die Verzerrungskomponenten &;: 


Eh +3 


und für die Änderung y,; des ursprünglich rechten Winkels zwischen den Koordinaten- 
flächen = const und «aj. = const: 


1 


Unter Berücksichtigung der Beziehungen: 
N 
E G Ir 
berechnet man leicht die g;, und daraus die Verzerrungen e, bzw. e, in der Richtung der 


Linien ®» =const bzw. u = const sowie den Winkel w, den die «- und v-Linien der verzerrten 
Mittelfläche einschließen: 


£ 


1 


wN, 


CO8@=SIN Ya Ya 


Dann ist 


11. Zur Ermittlung der Fundamentalgrößen L, M, N der deformierten Mittelfläche be- 
dienen wir uns erneut des absoluten Parallelismus und gelangen damit zu einer einfacheren 
und durchsichtigeren Herleitung wie die etwa bei Love mit Hilfe beweglicher Koordinaten- 
systeme gegebene'”). Dabei stehen zur Berechnung dieser Größen noch zwei Wege zur 
Verfügung: 

1. Zunächst können nochmals die Gleichgewichtsbedingungen nach dem in Nr. 6 und 7 
gegebenen Verfahren aufgestellt werden, indem man unter Berücksichtigung der Annahmen 
über den Deformationszustand der Mittelfläche (Nr. 10 und 13) von der Darstellung (15) der 
Mittelfläche ausgeht. In die so erhaltenen 6 Gleichungen des Gleichgewichts gehen jetzt statt 
der Größen Z, M, N — wie in den GIn. (13) und (14) — Ausdrücke in den Größen L, N 
und den Verschiebungskomponenten , r, w ein. Durch Vergleich entsprechender Glieder 
lassen sich Z, M, N bestimmen. 

2. Wir wollen jedoch hier einen anderen noch einfacheren Weg einschlagen, zu dessen 
Durchführung wir folgenden Hilfssatz benötigen: 

Der Vektor 

wird einmal vom Punkte P (u, ») der (deformierten oder undeformierten) Schalenmittelfläche 
im euklidischen Raume nach dem Nachbarpunkt P, (wu + du,» + de) der Schalenmittelfläche 
parallel verschoben. Dann wird derselbe Vektor in der in Nr. 6 konstruierten M, auf der 
Mittelfläche der Schale nach demselben Nachbarpunkt i im absoluten Sinne parallel verschoben. 
Die Komponenten der beiden so erhaltenen Vektoren in P, unterscheiden sich nur um Größen 
zweiter und höherer Ordnung in du, de. 

Man beweist diesen Satz leicht durch Rechnung, indem man die Gln. (12) auf die 
Komponenten des Vektors r anwendet. 


11) Siehe z. B. Trefftz: Mathematische Elastizitätstheorie in: Handbuch der Physik, Berlin 1928, Bd. VI, 8. 78. 
12), Love: |]. ec. 32%. 
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12. Nunmehr wird der Vektor v» im Punkte P(uw,v) der undeformierten Mittelfläche auf 
den zugeordneten Punkt P’(u,r) der «deformierten Mittelfläche übertragen, indem seine 
Komponenten auf das System X, X, X. (s. Gl. (15) und (19)) transformiert werden. Hierauf 
wird der transformierte Vektor vom Punkte P’(w,v) nach dem Nachbarpunkte P/ (u +du, v) 
der deformierten Mittelfläche parallel verschoben und man erhält den Vektor v/ in P/. In 
die Differentialgleichungen (12) der Parallelverschiebung gehen E, @, L, M, N ein. Dann 
wird derselbe Vektor » auf der undeformierten Mittelfläche vom Punkte P (u, v) zum Nachbar- 
punkte P,(u+du,v) parallel verschoben; in die Differentialgleichungen der Parallelver- 
schiebung gehen nur die bekannten Größen E, @, L, N ein. Schließlich wird der Vektor 

vom Punkte P, in den P, zugeordneten Punkt 

P!(u-+du,v) der deformierten Mittelfläche über- 

Pia») Pfurdu,v) tragen, indem seine Komponenten vom System 

| +du,e), (a +du,r), v(u+du,v) auf das 

Gl. 15) transformiert werden. Man erhält so den 

Pfiu) eformierte Mittelläche v, in Insgesamt wird also der Vektor », 

. EGIHN wie in Bild 5 symbolisch angegeben, auf zwei ver- 

Bild 5. schiedenen „Wegen“ I, II vom Punkte P (u, v) der un- 

deformierten Mittelfläche zum Punkte Pi(u+du,v) 

(ler deformierten Mittelfläche gebracht. Nach dem Hilfssatz in Nr. 11,2 unterscheiden sich 

die Komponenten (ler Vektoren v, und »; nur um Glieder von zweiter und höherer Ordnung 

in d«, da die Übertragungen P>P', P,>P‘; (Bild 5) nur den Charakter euklidischer Ko- 
ordinatentransformationen haben. 

Ist daher dt, der Zuwachs der ersten Komponente von » (auf das System X, X, Xu 
bezogen) nach Durchlaufung des „Weges“ I (Bild 5), d?, der entsprechende Zuwachs nach 
Durchlaufung des „Weges“ II, so folgt durch Grenzübergang d«—>0 mit lim (1)= en 


i 


au>o\du) 


du 


(Hiervon ist die partielle Ableitung t (u, zu unterscheiden!) 


Dabei geht in ga unter den Annahmen der Gln. (22) die unbekannte Größe 7, allein ein, im 


übrigen kommen wie auch in 3. nur die Fundamentalgrößen der undeformierten Mittelfläche 
und die Verschiebungen x, v, w vor, so daß /, errechnet werden kann. Insgesamt lassen sich 
so die Größen L, M, N bestimmen, wenn man dieselben Operationen auch für einen Punkt 


P,(u,e+de) und P}(u,ve-+do) der undeformierten bzw. deformierten Mittelfläche ausführt. 


13. Für die Übertragung des Vektors » vom x-X,-v-System im Punkte P (u, v) auf das 
Xu im Punkte P’ (w,v) gelten die Gleichungen: 


Die Verschiebungskomponenten «, ?, w müssen dem in Nr. 10 für die Deformation der Schale 
entworfenen physikalischen Modell entsprechen: 
1. Die Linienelemente senkrecht zur Mittelfläche sollen gerade bleiben und auch nach 
der Deformation senkrecht zur Mittelfläche sein. Das führt auf die Bedingungen: 


Verlangt man ferner, daß in der Richtung senkrecht zur Mittelfläche ki Dehnung oder 
Verkürzung erfolgt, so gilt weiter: 


Die folgenden Rechnungen sollen ferner zur Vereinfachung zunächst für den Fall einer 
dehnungslosen Deformation der Mittelfläche durchgeführt werden, so daß auch: 


(&,X,)=6, 2). 
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Doch lassen sich die ganzen Entwicklungen auch ohne die Bedingungen (22) durchführen 
(Nr. 14). 

Aus den Gln. (20) bis (22) folgen 6 Relationen zwischen den Verschiebungskomponenten 
x, r, w und ihren partiellen Ableitungen, so daß sich die Gln. (19) wie folgt vereinfachen: 


Vu 


al vN 

Überträgt man jetzt den Vektor v (P) auf die RGPOn Mittelfläche, so ergibt sich 


al 


Aus den Differentialgleichungen (12) der Parallelverschiebung folgt für die Verschiebung 
nach P!: 


E E. 


Daher wird mit lim Ii—°%. 
aher wird mı rc 


Andererseits ergibt die Parallelverschiebung von » (P) auf der undeformierten Mittelfläche 
(Weg II in Bild 5): 


v,(P,) 


in das System X, (u+du,r), X,(u+du,v), X„(u +du,v) übertragen: 


w u al al, 


Hieraus und aus Gl. (23): 


wu, 


dt, und dt, unterscheiden sich nach dem Hilfssatz Nr. 11,2 nur um Größen von der Ordnung 
du’, daher liefert jetzt der Grenzübergang d« >0 (Nr. 12): 


du du EYE (25). 


Setzt man (25) in (24) und daselbst für f,, ft, ihre Werte aus Gl. (23) ein, so wird endlich: 


Dann bestimmt man in gleicher Weise dt, und dt, und führt endlich dieselben Rech- 

nungen auch für die Punkte P,(u,v+dv) der undeformierten und P}(u,v+de) der defor- 
mierten Mittelfläche'?) aus. So erhält man: 


=, 1- Gu 1 Ev 


1m) Zur Rechnung werden jeweils die Komponentenänderungen dt, und dta benutzt, da in den Gleichungen für 
dt; immer 2 der unbekannten Größen Z, W, N gleichzeitig auftreten. 


(27), 


(28). 
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Die exakte mathematische Behandlung eines allgemeinen Problems der Schalentheorie 
stellt sieh nunmehr wie folgt dar: 


In den GIn. (13), (14) werden an Stelle der Spannungsresultanten und -momente die 
Verschiebungen eingeführt und die Fundamentalgrößen der deformierten Mittelfläche durch 
die Fundamentalgrößen der undeformierten Mittelfläche und die Verschiebungen gemäß den 
Gln. (19) und (26) bis (28) ausgedrückt. Aus dem so erhaltenen System partieller Differential. 
gleichungen sind die Verschiebungen unter Berücksichtigung der Randbedingungen zu be- 
stimmen. Dann berechnen sich die Dehnungen der Mittelfläche und alle die deformierte 
Mittelfläche geometrisch kennzeichnenden Größen aus den Gln. (19) und (26) bis (28). 


14. Werden die Rechnungen der Nr. 13 ohne die vereinfachende Annahme einer dehnungs- 
losen Deformation der Mittelfläche durchgeführt, so fallen die Gln. (22) weg. Die Gleichungen 
der Parallelverschiebung auf der deformierten Mittelfläche enthalten jetzt die metrischen 
Fundamentalgrößen E’, F’, @' (Nr. 10, Gl. 18). An Stelle einer unmittelbaren Bestimmung der 
zweiten Fundamentalgrößen L, M, N ergeben sich jetzt auf dem in den Nr. 12, 13 ein- 
geschlagenen Weg folgende 4 Relationen zwischen diesen Größen: 


+(@' EL + F' 


du’ 


3F'@ L+2E @ N+EF'N=(F'@ E@' +(F'@' Ei+ 


+2E6” (ü G du 


E )+2#® - 


F'@ L+2E @ G,+2E'G'F,— @' 6G,)t, 
N 

+2E”G 


of, 


E 

G)t +(E' 6,— 
dt, 

dv 


)+2R" @’- 


+2E 


w,, 
Aus irgend drei von diesen Gleichungen lassen sich L, M, N bestimmen. Für hin- 


reichend kleines F’, so daß die Produkte F’t; Größen zweiter Ordnung darstellen, sind die 
angenäherten Lösungen des Systems (29) gegeben dureh: 


Das im Abschnitt III angegebene Verfahren zur Berechnung der Größen Z, M, N läßt 
sich sogar noch anwenden, wenn auch die Bedingungsgleichungen (20), (21) wegfallen oder 
durch solche andere ersetzt werden, die sich aus einer entsprechenden Abänderung der An- 
nahmen über den Deformationszustand der Schale (Nr. 10 und 13) ergeben. Dieser Weg ist 
daher von größerer Allgemeinheit als die Berechnung der Verbiegungsgrößen mittels des 
Castiglianoschen Prinzips nach Trefftz'*) und führt mit geringerem Rechenaufwand zum 
Ziel. Andererseits liefert er aber den Zusammenhang zwischen den Spannungsresultanten 
bzw. -momenten und den Verschiebungen, die Trefftz gleichzeitig erhält, nicht mit, doch 
können diese natürlich auch in anderer Weise gewonnen werden. 340 


14) Trefftz: Z. angew. Math. Mech. Bd. 15 (1935), S. 108. 
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Spannungsverteilung in stabförmigen Membran- 
Kegelschalen. 


(Bericht der Festigkeitsabteilung der Focke-Wulf Flugzeugbau G.m.b.H.) 
Von A. Pflüger in Bremen. 

Bei einem biegungsfesten, geraden Stab, der einen dünnmwandigen Hohl- 
querschnitt mit längs der Stabachse veränderlichen Abmessungen besitzt, 
ist es für die Ermittlung der Spannungsverteilung nicht mehr ausreichend 
genau, wenn man den Stab an jeder gerade betrachteten Stelle durch einen 
Zylinder annähert und die hierfür gültigen Formeln der technischen Balken- 
biegungslehre benutzt. In der vorliegenden Arbeit wird die Spannungs- 
rechnung dadurch verfeinert, daß der Stab durch eine Kegelschale ersetzt 
wird. Die Ergebnisse der Theorie werden mit Messungen an einem Ver- 
suchsstück verglichen. 


1. Einleitung. 

Es ist allgemein üblich, bei der Spannungsberechnung biegungsfester, gerader Stäbe die 
Näherungstheorie der technischen Balkenbiegungslehre auch dann anzuwenden, wenn der 
Querschnitt des Stabes über die Stablänge veränderlich ist. Man ersetzt dabei den wirk- 
lichen Stab an jeder gerade betrachteten Stelle durch einen zylindrischen Stab gleichen Quer- 
schnitts. Dieses Verfahren ist jedoch dann nicht mehr zulässig, wenn der Stab einen dünn- 
wandigen Hohlquerschnitt besitzt. Hierbei zeigt sich nämlich, daß die auftretenden Schub- 
spannungen, die bei Hohlquerschnitten im allgemeinen sehr wesentlich sind, schon recht er- 
heblich beeinflußt werden, wenn die Form des Stabes auch nur wenig von der eines Hohl- 
zylinders abweicht. Diese Tatsache läßt sich für den bisher in der Literatur behandelten 
Sonderfall, daß der Stab aus einem Kreis-Hohlkegel besteht, leicht nachweisen). Der Gegen- 
stand der folgenden Untersuchungen soll die Spannungsverteilung sein, die in einem stab- 
förmigen Hohlkegel allgemeiner Art bei Beanspruchung durch Längskraft, Querkraft, Biegung 
und Torsion hervorgerufen wird. Um mit der gestellten Aufgabe eine anschauliche Vorstellung 
zu verbinden, denke man bei den folgenden Rechnungen etwa an den in Bild 1 dargestellten 
einseitig eingespannten Balken, der an seinem freien Ende durch Kräfte und Momente belastet 
ist. Wenn der Hohlkegel auch noch nicht die allgemeinste Form ist, die ein Stab mit Hohl- 
querschnitt haben kann, so wird man doch stets eine bessere Lösung erhalten, wenn man 
den zu berechnenden Stab durch einen Kegel annähert, als wenn man ihn nur durch einen 
Zylinder ersetzt. 

2. Allgemeines, Bezeichnungen. 

Bild 2 zeigt einen Teil der Mittelfläche der zu untersuchenden Kegelschale. Dabei wird 
als Stabachse eine Gerade bezeichnet, von der vorläufig nur vorausgesetzt sei, daß sie durch 
die Kegelspitze hindurchgeht und innerhalb der Schale verläuft; ihre genaue Lage wird erst 
später festgelegt. Von einem stabförmigen Kegel kann nur dann gesprochen werden, wenn 
seine Abmessungen in Richtung der Stabachse wesentlich größer sind als in den beiden dazu 
senkrechten Richtungen. Der Winkel zwischen der Kegelerzeugenden und der Stabachse muß 
also stets ein kleiner Winkel sein, der die Größenordnung von 10° nicht allzusehr über- 
schreiten darf. Die Schnittkurve, die die Ebene eines Querschnitts senkrecht zur Stabachse 
mit der Schalenmittelfläche bildet, sei als Querschnittskurve bezeichnet. Sie wird im all- 
gemeinen keine Krüm- 
mungslinie der Mittel- 
fläche sein, jedoch ‚we- 


Drehmoment gen der stabförmigen 
fer Gestalt des Kegels nicht 
öngskraft 


sehr davon abweichen. 
Der Einfachheit halber 
ist in den Bildern 1 und 2 
ein Kegel mit „einzelli- 
gem“ Querschnitt dar- 
gestellt worden. Die 
folgenden Betrachtun- 
gen gelten jedoch auch 
für mehrzellige Quer- 
schnitte. 


Querkraft 


Bild I (oben). Einseitig eingespaunter 


kegelförmiger Stab mit Hohlquersehnitt. Krümmungslinie 


Bild 2 (rechts). Bezeichnungen bei der 
Kegelschale. Querschnittkurve 


HM. Neuber: Z. angew. Math. Mech. Bd. 14 (1934), S. 203. — H. Ebner: Luftf.-Forschg. Bd. 14 (1937), S. 9%. 
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Es mögen bedeuten 


x, 9, 2 Koordinaten eines rechtwinkligen Systems, dessen Nullpuukt in der Kegelspitze 

liegt, und dessen x-Achse mit der Stabachse zusammenfällt, 

s Abstand eines Punktes der Schalenmittelfläche von der Kegelspitze, 

a Winkel zwischen der Stabachse und den Erzeugenden der Mittelfläche, 

o Krümmungsradius der Querschnittskurven, 

y in der Querschnittsebene gemessener Winkel, den ein Krümmungsradius o mit 
einer (beliebig wählbaren) Nullrichtung bildet, 

© Winkel zwischen den Querschnittskurven und den Krümmungslinien der Schalen- 
mittelfläche, die zu den Erzeugenden des Kegels orthogonal sind. 


Die als positiv gewählte Richtung von o, 9 und & geht aus Bild 2 hervor. Die Gestalt 
‚der Mittelfläche ist festgelegt, wenn für einen bestimmten Querschnitt & —=konst. der Verlauf 
der Querschnittskurve gegeben ist. Im folgenden seien alle Abmessungen der Schale auf 
einen Schnitt hezogen, der sich im Abstand «=. von der Kegelspitze befindet. Es sei 


Ya; 2a Koordinaten der Querschnittskurve des Schnittes a, 
0a Krümmungsradius der Querschnittskurve des Schnittes z= a. 


3. Gleichgewicht am Stabelement. 


In der Theorie der gewöhnlichen Balkenbiegung werden nur solche Spannungszustände 
betrachtet, die sich in hinreichender Entfernung von einer Lastangriffsstelle einstellen. Die 
durch die Art der Krafteinleitung bedingten Spannungsstörungen sind nach dem Saint- 
Venantschen Prinzip örtlich beschränkt und außerdem im allgemeinen praktisch von unter- 
geordneter Bedeutung. Bei dem hier zu berechnenden Kegelstab soll die Spannungsverteilung 

unter sinngemäßer Erweiterung der bekannten Ansätze der 
Balkenbiegung ermittelt werden. Es sei daher zur Unter- 
Bild 3. Schnittkräfte suchung des Gleichgewichts an einem Stabelement, das durch 
zwei Schnitte e=konst. aus dem Kegel herausgeschnitten 
j + wird, ebenfalls eine Stelle des Stabes betrachtet, an der 
keine äußere Belastung wirkt. An dem in Bild 3 dargestellten 
Stabelement greifen demnach nur Schnittkräfte bzw. -mo- 
mente an. Es sei 


P Längskraft, 

Querkraft parallel zur z-Achse, 

ce Abstand der Querkraft Q von der «-z-Ebene (über 
die Größe von e wird später noch besonders verfügt), 

M Biegemoment in der x#-z-Ebene (Momentenvektor 
parallel zur y-Achse), 

D Drehmoment. 


Die positive Richtung dieser Größen geht aus Bild 3 

hervor. Zur Vereinfachung der folgenden Rechnungen soll 

eine Beanspruchung durch Querkräfte und Biegemomente in der x-y-Ebene nicht berück- 
sichtigt werden. Sie kann unbeachtet bleiben, da ihre Untersuchung nichts Neues bietet und 
eine spätere Verallgemeinerung der gewonnenen Ergebnisse keinerlei Schwierigkeiten bereitet. 


Die Bedingungen für das Kräftegleichgewicht in Richtung der x- und z2-Achse und für 
das Momentengleichgewicht um die x-Achse liefern das in Bild 3 bereits vorweg genommene 
Ergebnis, daß die Längskraft P, die Querkraft 0 und das Drehmoment D unabhängig von x 
einen konstanten Wert haben müssen. Aus der Momentenbedingung um eine zur y-Achse 
parallele Gerade folgt 


M=M,+ 0x, 


wobei die Integrationskonstante M, als ein in der Kegelspitze angreifendes Moment gedeutet 
werden kann. M, sei durch die Forderung festgelegt, daß das Biegemoment an der Stelle 
x—.a den Wert M= M,„ annehmen soll. Damit wird 


Muh? 
y 
M+dM 
_ 
XUM 
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4. Gleichgewicht am Schalenelement. 


Während Gl. (1) etwas über das Gleichgewicht „im Großen“ aussagt, sollen jetzt die 
Verhältnisse am einzelnen Schalenelement betrachtet werden. Dabei sei die wesentliche ver- 
einfachende Annahme gemacht, daß die Schalendicke hinreichend gering ist, um den Kegel 
als Membranschale, d. h. als Schale ohne Biegungs- und Drillungssteifigkeit der Mittel- 
fläche, berechnen zu können. Wenn es sich auch im folgenden um das Aufsuchen eines 
Spannungszustandes handeln wird, der wie die Ansätze der gewöhnlichen Balkenbiegung im 
allgemeinen nur eine Näherung darstellt, so soll doch von der gesuchten Spannungsverteilung 
verlangt werden, daß sie die Gleichgewichtsbedingungen am Element der Kegelschale exakt 
erfüllt. 


Zur Aufstellung dieser Bedingungen sei aus der Schalenmittelfläche ein Element durch 
zwei benachbarte Schnitte «=konst. und 9=konst. herausgeschnitten. Die Begrenzungs- 
linien des Elementes sind also zwei Querschnittskurven und zwei Erzeugende; die Seiten- 
längen sind od und ds. Da die Querschnittskurven und die Erzeugenden auf der Schalen- 
mittelfläche ein Netz schiefwinkliger Koordinatenlinien darstellen, werden an dem betrachteten 
Schalenelement auch am zweckmäßigsten „schiefe Schnittkräfte“ definiert ?). 


In Bild 4 ist das Schalenelement mit den daran 
angreifenden Kräften dargestellt. Die in Richtung der Er- En 
zeugenden auftretenden Längskräfte sind mit N, die Schub- u Een 
kräfte (oder Schubflüsse) mit 7 bezeichnet. Eine äußere j 
Belastung ist aus den schon erwähnten Gründen auch 
hier nicht vorhanden. Querkräfte, Biegemomente und 
Drillungsmomente treten nicht auf, da es sich voraus- 
setzungsgemäß um einen Membranspannungszustand han- 
delt. Schließlich ist noch zu erwähnen, daß in Bild 4 
eine in Richtung der Querschnittskurve wirkende Längs- 
kraft nicht eingetragen ist. Ihr Vorhandensein würde an 
sich mit der Forderung eines Membranspannungszustandes 
verträglich sein. Aus der in Mgdy)as 
Richtung der Normalen zur Schalenmittelfläche folgt je- 
doch sofort, daß die in Frage stehende Schnittkraft im 
vorliegenden Fall verschwinden muß. Sie würde nämlich 
allein einen Beitrag zu dieser Gleichgewichtsbedingung liefern, da die Kräfte N und 7 in der 
Tangentialebene liegen und eine äußere Belastung des Schalenelementes nicht vorhanden ist. 
Von den zur Verfügung stehenden Gleichgewichtsbedingungen ist damit eine bereits ver- 
braucht. Die Momenten-Gleichgewichtsbedingung um die Schalennormale sagt ebenfalls nichts 
Neues mehr aus; sie liefert die in Bild 4 schon als bekannt vorausgesetzte Gleichheit der 
einander zugeordneten Schubkräfte. Zur Aufstellung von Beziehungen zwischen den Schnitt- 
kräften N und 7 bleiben dann noch zwei Bedingungen für das Kräftegleichgewicht in der 
Tangentialebene übrig. Zweckmäßig wird hier das Gleichgewicht in Richtung der Tangente 
an die Krümmungslinie und in Richtung der Erzeugenden betrachtet. Als erste Bedingung 
erhält man dann 


(To dp)dscos»+ Tds 


dabei ist cos do der Winkel, den die beiden das Schalenelement begrenzenden Erzeugen- 


den miteinander bilden. Als Gleichgewichtsbedingung in Richtung der Erzeugenden ergibt sich 


y)dssino=®. 


Beachtet man, daß cos«a -- ist, und setzt zur Abkürzung 


d(...) 


so ergibt sich aus der ersten der beiden Gleichgewichtsbedingungen nach Division durch dedo 


2) Flügge: Statik und Dynamik der Schalen, Berlin 1934, S. 11. 
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und aus der zweiten 


cos a 


(No) + (To) sin®—=0 
oder 


Es ist nun 


Damit wird 


Bezieht man wieder alle Ausdrücke auf den Querschnitt z-—=a, indem man 0=x 
setzt, so nehmen die beiden Gleichgewichtsbedingungen, denen der gesuchte Spannungszustand 
genügen muß, folgende Form an 
(Ta) +T 
2a, b). 
cosa+ T=0 (2a, b) 
Sind die Schnittkräfte N und T bekannt, so ist es z. B. für die Berechnung der in der 
Schale auftretenden Spannungen erforderlich, die schiefen Schnittkräfte auf die Schnittkräfte 
eines rechtwinkligen Netzes von Koordinatenlinien der Schalenmittelfläche zu transformieren. 
Hıerzu wird zweckmäßig das System der Krümmungslinien gewählt. Die hierbei in Richtung 
der Erzeugenden auftretende Längskraft sei N* und die zu- 
gehörigen Schubkräfte 7T*. Die Längskraft in Richtung der 
zur Erzeugenden orthogonalen Krümmungslinie muß wieder zu 
Null werden, da sonst die Gleichgewichtsbedingung in Richtung 
ee der Schalennormalen nicht erfüllt werden kann. Zur Aufstellung 


—ds 


ı 177) der Beziehungen zwischen N, Tund N*, T* wird nach Bild 5 ein 
ds | Querschnitkurre  dreieckiges Element der Schalenmittelfläche betrachtet, dessen 

Tds TE Seiten der Erzeugenden, der Tangente an die Querschnittskurve 
Erzeugende N und der Tangente an die zur Erzeugenden senkrechte Krüm- 
mungslinie parallel sind. Aus der Gleichgewichtsbedingung in 


Bild 5. Dreieckiges Schalenelement. s-Richtung und senkrecht dazu folgt dann sofort 


(N + Tsin ®) +Tds NM*dsetgo 


Tcos® — — 
sin 


Oder aufgelöst nach N* und T* 


N 
N'= 2T tg 


"=? 


5. Ansatz für die Spannungsverteilung. 


‘ Durch die bisher aufgestellten Formeln’ ist die Spannungsverteilung durchaus noch nicht 
festgelegt. Vielmehr lassen sich beliebig viele Spannungszustände angeben, die sämtlich die 
bisher allein betrachteten Gleichgewichtsbedingungen erfüllen. Versucht man zu einer 
Festlegung des in Frage kommenden Spannungsansatzes dadurch zu gelangen, daß man auch 
die Verformungsbedingungen der Schale betrachtet, so ergibt sich zunächst folgende 
Schwierigkeit, die nicht nur für die Kegelschale,- sondern auch schon für eine zylindrische 
Schale mit dünner Wandstärke besteht. Es war vorausgesetzt worden, daß der zu unter- 
suchende Stab als Membranschale berechnet werden kann. Bei Membranschalen führt aber 
im allgemeinen jeder die Gleichgewichtsbedingungen befriedigende Spannungszustand auch zu 
einem möglichen Verformungszustand. Die Festlegung der sich bei einer bestimmten Be- 
lastung einstellenden inneren Kräfte der Schale geschieht erst durch die besonderen Rand- 
bedingungen, die für den Angriff der äußeren Kräfte vorgeschrieben sind. Der in dem Stab 
auftretende Spannungszustaud würde damit von der Art der Einleitung der Kräfte P und @ 


(N o)' r sino 
| 
= ( 1 ) 
o \cos a 
XUM 
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1942 


und der Momente M und D auch an den Stellen abhängig werden, die vom Ort der Ein- | 
leitung weit entfernt sind®). Es würde dann die ganze beabsichtigte Übertragung der Ge- 
dankengänge der gewöhnlichen Biegelehre auf Stäbe mit dünnwandigem Hohlquerschnitt über- 
haupt nicht mehr möglich sein. Die Verhältnisse liegen jedoch in Wirklichkeit dadurch anders, 
daß eine reine Membranschale, wie sie für die bisherige Untersuchung vorausgesetzt war, 
praktisch nie vorliegt, da sie wegen ihrer Fähigkeit, große spannungslose Verformungen aus- 
führen zu können, gar kein brauchbares Bauelement sein würde. Vielmehr ist es stets er- 
forderlich, der Schale durch Rippen oder Spante und im allgemeinen auch durch Längs- | 
versteifungen in Richtung der Erzeugenden eine gewisse Steifigkeit zu geben, falls hierzu | 
nicht schon die bisher vernachlässigte, in Wirklichkeit aber stets vorhandene Biegesteifigkeit 
der Schalenmittelfläche ausreichend ist. Infolgedessen kann auch der Stab mit dünnwandigem 
Hohlquerschnitt nur noch ganz bestimmte Verformungen ausführen, die denen eines Stabes 
mit Vollquerschnitt genau entsprechen. Die an Krafteinleitungsstellen auftretenden Spannungs- 
störungen werden durch die Versteifungen bzw. durch den Biegewiderstand der Schalen- 
mittelfläche zum Abklingen gebracht. Zur Aufstellung eines Spannungsansatzes können dann 


nd die bekannten Annahmen der Balkenbiegung, z. B. die Annahme vom Ebenbleiben der Quer- 
schnitte, auch bei stabförmigen Membranschalen sinngemäß angewendet werden. Der so ge- 
wonnene Spannungsansatz wird jedoch — das ist das wesentliche Ergebnis der vorstehenden 
Überlegungen — ebenso wie beim Stab mit Vollquerschnitt im allgemeinen nur eine Näherung 
b). sein, da die zu dem angenommenen Membranspannungszustand exakt gehörenden Formände- | 
rungen wegen der vorhandenen Versteifungen oder wegen der Biegesteifigkeit der Schale 
nicht möglich sind. 
ler Bei der Kegelschale sei mit dem Aufsuchen eines Ansatzes für die Schnittkräfte N be- 
fte gonnen. Der Verformungszustand sei durch die Annahme festgelegt, daß bei torsionsfreier 
en; Beanspruchung aus Längskraft, Querkraft und Biegung ein zur Stabachse senkrechter Quer- 
ng schnitt bei der Deformation eben bleibt. Die auf Grund dieser Annalıme gewonnene Lösung 
zo wird bekanntlich bei zylindrischen Stäben mit der exakten Lösung identisch, wenn der Stab 
er eine reine Biegungsbeanspruchung erfährt. Für kegelförmige Stäbe wird dieses jedoch nicht 
m mehr gelten. Für den einfachen Fall reiner Längskraftbeanspruchung durch eine Zugkraft P 
ng ist in Bild 6 ein durch zwei Schnitte « = konst. begrenztes Stabelement vor und nach der 
ng Verformung dargestellt. Die Zeichenebene ist dabei eine die Kegelachse enthaltende Schnitt- 
ebene. Es sei 
en P 
ve ö Wandstärke der Schale, 
m- E Elastizitätsmodul, 
ın &; Dehnung der Schale in Richtung der 
Erzeugenden, 
Änderung des Abstandes dx der beiden 
benachbarten @uerschnitte des Stab- 
A Bild 6. Stabelement der Kegelschale bei Beanspruchung 
elementes. dureh Längskraft. 
Nach Bild 6 ist dann 
Bei Vernachlässigung quadratischer Glieder der Dehnungsgrößen wird 
| = 
ds 
Da in der Schalenfläche keine Längskräfte in Richtung der Krümmungslinien wirken, nimmt 
ht das Elastizitätsgesetz für die Schnittkräfte N hier die einfache Form an 
ie N*=Ede,. 
2 Nach (3a) ist bei Abwesenheit von Schubkräften 
le N 
he TFT 
T- 
er Man erhält also für die Schnittkräfte X 
N=AEdcosweos’a. 
— 
d- 3) Diese Tatsache widerspricht keineswegs dem Saint-Venantschen Prinzip, das nur dann gilt, wenn die 
| Gleichgewiehtsgruppen, über deren Abklingen etwas ausgesagt wird, in einem hinreichend kleinen Gebiet angreifen. 
ıb Die gesamte Länge der Querschnittskurve, über die sich die hier betrachteten Spannungen verteilen, ist eben schon ein 


zu großer Bereich der Schalenfläche. 
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Das gilt nur für eine Beanspruchung durch eine Längskraft P. Wirkt außerdem auf 
das Stabelement ein Biegemoment M, so ist der obige Ansatz für N der Annahme vom Eben- 
bleiben der Querschnitte entsprechend zu erweitern: 


B ist dabei wie A eine Konstante des betrachteten Querschnitts. Würde der Stab auch durch 
ein Biegemoment in der x-y-Ebene beansprucht werden, so müßte in (4) in der Klammer noch 
ein weiteres Glied hinzugefügt werden. Der Einfluß von Verformungen, die durch Schub- 
kräfte hervorgerufen werden, wird natürlich durch den Ansatz (4), wie üblich, vernachlässigt. 
Selbstverständlich bleibt noch zu prüfen, ob trotz aller Annahmen und Vernachlässigungen 
die Beziehung (4) zu einer Spannungsverteilung führt, die die Gleichgewichtsbedingungen (2) 
widerspruchsfrei erfüllt. 


Zur Bestimmung der Konstanten A und B dienen die Bedingungen, daß die im Quer- 
schnitt übertragenen Schnittkräfte eine Resultierende von der Größe P und ein resultierendes 
Moment von der Größe M liefern müssen: 


PNeosaody = 
PNeosazody=M. 


N cos a ist dabei die auf der Querschnittsebene senkrecht stehende Komponente der 
Kräfte N. Das Zeichen p soll andeuten, daß das Integral über die ganze Querschnittskurve, 
d. h. über alle Winkel 9 zu erstrecken ist. Im folgenden sei 


F=9dody, J=9öz’ody, 
ö=6c0s w cos’a. 
F=9Ödodp, J=$Öözodp. 


F und J ergeben sich als Fläche und Trägheitsmoment des Querschnittes an der Stelle «, 
wenn jedem Punkte der Querschnittskurve die Wandstärke ö zugeschrieben wird. Nach Be- 
rechnung von A und B wird dann aus (4) 

M 
a 


F + 2 
unter der Voraussetzung, daß für die Lage der Achsen y und z die bekannten Bedingungen 
Pöyodp=0, Pözody=0, Pöyzoedy=0 


vorgeschrieben werden. Da diese Beziehungen für jede Stelle x gelten sollen, wird dadurch 
die Stabachse als Verbindungslinie der Schwerpunkte der mit der Wandstärke ö belegten 
Querschnittskurven definiert. Ferner ist jetzt die Lage der Achsen y und z festgelegt, die 
den Hauptachsen des Querschnitts mit der Wandstärke ö parallel sein müssen. Schließlich 
ist noch durch die obigen Bedingungen über den Verlauf der Wandstärke eine bestimmte 
Voraussetzung getroffen. Damit nämlich die Querschnittsschwerpunkte auch tatsächlich auf 
einer durch die Kegelspitze gehenden Geraden liegen, muß 


fı(p) 


sein, wobei f eine reine Funktion von x und f, eine reine Funktion von p ist. Für den 
Querschnitt an der Stelle = a sei 


(5) 


d=6,(P), fl, F=F,, F 
Es gilt nun 


I= 
so daß man aus (5) unter Benutzung von (1) erhält 
da 


P 


Die Größe der Längskräfte N ist damit für jede Stelle der Schale festgelegt. 


| 
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6. Die Schubkräfte. 


Es sei jetzt nachgewiesen, daß mit dem gewählten Ansatz (5) die Gleichgewichts- 
bedingungen (2) erfüllt werden können. Hand in Hand damit geht die Bestimmung der 
Schubkräfte T. Setzt man (5a) in (2b) ein, so folgt 


a? 
Zur Abkürzung sei gesetzt 


Man erhält dann 


a? 


T=C(&) — aZa0adY, 


wobei € eine nur von x abhängige Funktion ist. Setzt man den Ausdruck für T in (2a) ein, 
so wird 
+0 = 
a 
(= Co x? 


wenn man vorschreibt, daß C für »—=a den Wert C, annehmen soll. Man erhält dann für 
dıe Verteilung der Schubkräfte 


oder in der bequemeren Schreibweise, bei der jedoch die Abhängigkeit von x nicht mehr un- 
mittelbar zum Ausdruck kommt: 


0 
Artet der Kegel in einen Zylinder aus, so ergibt sich für den Schubfluß die bekannte 
Formel 
0 
die sich von (6) nur durch den fehlenden Querstrich über den Größen @, J und ö unter- 
scheidet. Da die Winkel a und » klein sind, werden J und Ö immer nur wenig von J und ö 


A 
abweichen. Der Betrag =. um den beim Kegel die Querkraft Q abgemindert werden muß, 


um Q zu erhalten, ist jedoch meist sehr wesentlich. Diese Abminderung erklärt sich an- 
schaulich daraus, daß die Schnittkräfte N Komponenten Nsina in der Querschnittsebene 


haben, die zusammen in Richtung von @ gerade den Betrag = ausmachen. Es ist nämlich 


2 
0 — N Sina 77. 5049, 
wobei vr E zZ der Kosinus des Winkels ist, den die Komponenten N sin a mit der 2-Achse 
bilden. Mit 


h. 
ıf 
). 
h 
)- . 
t. 
n 
) 
| 

) 

n 
h 
n 
e 
h 
e 
| 
sina= cos a 
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wird dann unter Benutzung von (5) 


M, 2 


was zu beweisen war. Der Angriffspunkt des Querkraftanteils z ist der Punkt #=0, y=V. 


Dieses ergibt sich sofort aus der Tatsache, daß alle Längskräfte N die Stabachse schneiden. 
Die Berechnung der Integrationskonstanten (', erfordert besondere Überlegungen. Die 


Größe C=(, = bedeutet für einen bestimmten Querschnitt einen konstanten Schubfluß, der 


weder in y- noch in z-Richtung eine resultierende Kraft, sondern lediglich ein resultierendes 
Drehmoment liefert. Gl. (6) gilt in der angegebenen Form nur für eine zweifach zusammen- 
hängende Querschnittsfläche, d. h. also für einen Stab mit einzelligem Querschnitt. Besitzt 
der Stabquerschnitt mehrere Zellen, so ist es mit den Gleichgewichtsbedingungen verträglich, 
wenn für jede Zelle eine besondere Konstante C als Schubfluß an der Stelle 9=0 vor- 
geschrieben wird. Die zur Ermittlung von € erforderlichen Rechnungen sollen jedoch auch 
weiterhin am einzelligen Querschnitt erläutert werden. Über die Verallgemeinerung der auf- 
gestellten Formeln braucht dann später nur noch wenig gesagt zu werden, da die hierzu not- 
wendigen Überlegungen die gleichen sind, wie in dem in der Literatur bereits mehrfach be- 
handelten einfacheren Fall des zylindrischen Stabes. Die Größe € wird zweckmäßig zerlegt in 


Die in einem Querschnitt übertragenen Schubkräfte müssen in bezug auf die Stabachse 
ein resultierendes Moment von der Größe D+ Ye ergeben. Es sei nun verlangt, daß die 
durch CP dargestellte Schubflußverteilung gerade das Drehmoment D liefert, während (9 zu- 


sammen mit den sonst noch vorhandenen von Q abhängigen Schubkräften das Moment ®e 
liefern soll. 


7. Torsionsfreie Verformung. Querkraftmittelpunkt. 


Über der Abstand e ist bisher noch nicht verfügt worden. Bei einer bestimmten ge- 
gebenen Belastung des Stabes ist D oder c frei wählbar. Im folgenden möge die Frage unter- 
sucht werden, wie groß (@ gewählt werden muß, damit e gerade die y-Koordinate des „Quer- 
kraftmittelpunktes“ ist. Bei zylindrischen Stäben wird dieser Punkt auch sehr häufig „Schub- 
mittelpunkt“ genannt. Hier soll jedoch der ersten Bezeichnung der Vorzug gegeben werden, 
um damit anzudeuten, daß es sich um einen Punkt handeln soll, in dem die gesamte im Quer- 
schnitt übertragene Querkraft Q und nicht etwa nur der von den Schubkräften T übernommene 
Beitrag Q angreift. Der Querkraftmittelpunkt ist derjenige Punkt, durch den @ hindurchgehen 
muß, damit der Stab an der betrachteten Stelle eine reine Biegungsverformung und keine 
Torsionsverformung erleidet. Diese Definition ist jedoch zur formelmäßigen Aufstellung der 
für die Berechnung von C@ notwendigen Bedingungen noch nicht ausreichend, da ja die 
Biegungsverformung im Abschnitt 5 lediglich angenähert charakterisiert worden ist und über 
die Torsionsverformung überhaupt noch nichts gesagt worden ist. Eine präzise Vorschrift über 
die Ermittlung von C9 erhält man, wenn man von den durch die Gl. (5) und (6) gegebenen 
Spannungszuständen ausgeht und zur Festlegung des Querkraftmittelpunktes nach Trefftzt) 
die Formänderungsarbeit benutzt°). Der zur reinen Torsionsverformung gehörende Spannungs- 
zustand sei durch den Index D, der zur reinen Querkraftbiegung gehörende durch den Index U 
gekennzeichnet. Es ist 


(CP; 


Es soll nun von torsionsfreier Biegung dann gesprochen werden, wenn die beiden 
Spannungszustände mit den Indizes D und @ zueinander orthogonal sind in dem Sinne, daß 
der eine Zustand auf den Verschiebungswegen des anderen keine Arbeit leistet. Bei gleich- 
zeitiger Beanspruchung durch Querkraft-Biegung und Torsion kann also die gesamte Form- 
änderungsarbeit additiv aus den Arbeiten der einzelnen Zustände zusammengesetzt werden. 


4) E. Trefftz: Z. angew. Math. Mech. Bd. 15 (1935), S. 222 


5) Ein anderer Weg wird von W. T. Koiter: Versl. en Verh, Nat. Luchtvaartlaborat. Amsterdam, Bd. 8 (1939), 
S. 169, eingeschlagen. 


ND-0, 
79= — 7 \dzedg. 
0 
XUM 
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Es sei zunächst allgemein die Formänderungsarbeit ermittelt, die ein aus den Sehnitt- 
kräften N, TI bzw. T*! bestehender Spannungszustand auf den Verschiebungswegen 
eines Zustandes NH, TU bzw. N*l, T*U je Einheit der Schalenmittelfläche leistet. Diese Arbeit 
sei mit u; bezeichnet und ist, wenn @ der Schubmodul ist, 


Go 


Weitere mit der Querkontraktion behaftete Glieder kommen nicht vor, da in Richtung der 
Krümmungslinien keine Längskräfte wirken. Nach Einführung der Schnittkräfte des schief- 
winkligen Schalenelementes ergibt sich 


1,11 1 N! 

dx, EÖ 
Ed Gö +4 o)+ Eö ' 


yıqu 


+27 + 


— 2 T!tg - 
cos © + \cos m 


Das schiefwinklige Schalenelement hat die Fläche odpdscosw=o dpdx. Die 


im folgenden mit 4%" bezeichnete Arbeit des Stabes je Längeneinheit der Stabachse an der 
Stelle x ergibt sich dann durch Integration über alle p zu 


1,1 _c0S@ 
x 90 


NINU 1 G TINU tego ]|coso 
4x $| Ed | E ig o)+ Eö cos cosa ® 


Setzt man in (7) 
NI=ND TI=TD, NU=N® Ti-70 
und verlangt, daß A’ zu Null wird, so folgt 


G 2 ö M 
0 


cosa” 


oder 


1 

3 | Özol +4 tg? o) cos’ cosal odg=d. 
0 - 

Aufgelöst nach ©® ergibt sich 


mit 


) G 


Fi +4 tg? o) cos’ wcos’aodp 


zsin 


Wird wieder alles auf den Querschnitt an der Stelle #—=« bezogen, so erhält man 


2 


a a 


k,=k,a 


1 
und damit 


a? 
c9= Vakıat [M. a)] 
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Zur anschaulichen Deutung der durch die Gl. (8) 
und (8a) gewonnenen Aussage wird zweckmäßig die Querknaftmitelgunkt 
Koordinate ce des Querkraftmittelpunktes berechnet. 

Bild 7 zeigt den Stabquerschnitt an der Stelle x. Längs 

der Querschnittskurve greifen die Schubkräfte 7TQ an. 

An Stelle der in der Querschnittsebene wirken- 

den Komponenten der Längskräfte N ist deren Re- vu 


v 


sultierende = eingetragen, von der bereits bekannt ist, 


daß sie durch die Stabachse hindurchgeht. Die im Abstand c angreifende Querkraft Q muß 
nun um die Stabachse dasselbe Moment liefern wie die Schubkräfte 7T®. Ist r das Lot vom 
Punkte 2<=0, y=0 auf die Tangente an die Querschnittskurve, so ergibt sich 


T9rodg 
—=(9k,+Mk)2% 


wenn 3-4 $ rodg die von der Querschnittskurve umschlossene Fläche bedeutet. Mit 


erhält man 


0 
und mit du 


Als wichtigste Aussage folgt aus (9a), daß c nicht wie bei zylindrischen Stäben eine 


Konstante ist, sondern eine von x und von dem durch die Belastung bedingten Verhältnis S 


linear abhängige Größe. Die Abhängigkeit von « bedeutet, daß nur für einen bestimmten 
Querschnitt und nicht für den ganzen Stab torsionsfreie Biegungsverformung erzielt werden 
kann. Da sich für die 2-Koordinate des Querkraftmittelpunktes, Gl. (9a) entsprechend, eben- 
falls eine lineare Abhängigkeit von x ergibt, ist die Verbindungslinie der Querkraftmittel- 
punkte, die Querkraftmittellinie, eine Gerade. Die Winkel, die diese Gerade mit den Rich- 
tungen der Koordinatenachsen bildet, sind bei gegebenen Abmessungen des Stabes konstant. 
Im übrigen ist jedoch die Lage der Geraden eine Funktion der Belastung. Besteht die Be- 
lastung des Stabes aus einer in der Kegelspitze angreifenden Einzellast, so wird @=0; die 
Querkraftmittellinie ist dann eine durch die Kegelspitze hindurchgehende Gerade. Wird je- 
doch 9=0, so geht e nach unendlich. Dieses erklärt sich durch die Tatsache, daß die quer- 
kraftfreie Biegung nicht torsionsfrei ist. Um bei Q=0 verdrehungsfreie Verformung zu er- 
zielen, ist es vielmehr notwendig, noch zusätzlich ein Drehmoment aufzubringen; hierzu ist 
bei im Grenzfall verschwindender Querkraft ein unendlich großer Abstand ec erforderlich. 
An die der Beziehung (8) zu Grunde liegenden Arbeitsgleichung sei noch folgende Be- 
merkung geknüpft. Man gelangt zu derselben Gleichung, wenn man vorschreibt, daß der 
konstante Schubfluß (2 gerade so groß gewählt werden soll, daß er auf den Verschiebungs- 
wegen des gesuchten Spannungszustandes N®, TQ keine Arbeit leistet. Wegen dieser Ortho- 
gonalitätsbedingung stimmt die durchzuführende Rechnung rein formal mit der Berechnung 
eines statisch unbestimmten Systems genau überein. Die Größe C®9 entspricht der Un- 
bekannten, die Schnittkräfte mit 09 =0 entsprechen der Spannungsverteilung am statisch 
bestimmten Hauptsystem. Die Analogie läßt sich an einem geeigneten Modell°) soweit ver- 
vollständigen, daß man dort auch die Verschiebungswerte, die bei der Berechnung des statisch 
unbestimmten Systems auftreten, anschaulich deuten kann. Diese Methode ist vor allem bei 
mehrzelligen Querschnitten zylindrischer Stäbe zweckmäßig; ihre Anwendung dürfte sich 
jedoch bei den komplizierteren Beziehungen, die für kegelförmige Stäbe gelten, weniger lohnen. 
Es sei jetzt kurz auf die Verhältnisse bei mehrzelligen Querschnitten eingegangen. Wie 
bereits erwähnt, sind die in diesem Falle notwendigen besonderen Rechnungen im Prinzip bei 
%) Vom Verfasser ist in Ing.-Arch. Bd. 8 (1937), S. 25, für zylindrische Stäbe vorgeschlagen worden, den Stabaquer- 


schnitt als ein Rahmensystem aufzufassen, das senkrecht zu seiner Ebene durch eine gedachte Belastung geeignet bean- 
sprucht wird. 
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einem Kegel dieselben wie bei einem zylindrischen Stab’). Die Längskräfte NQ des Spannungs- 
zustandes reiner Querkraft-Biegung sind nach wie vor durch Gl. (5) mit P=0 gegeben. Die 


Verteilung der Schubkräfte 79 ändert sich insofern, als jetzt in jeder Zelle eine Größe C9 an 
der Stelle unbekannt ist, wo der Weg des Integrals |öz0odp beginnt. Beim Spannungs- 


0 
zustand reiner Torsion wirkt in jeder Zelle ein für die einzelnen Zellen verschieden großer 
konstanter Schubfluß C?. Das resultierende Moment aller Schubkräfte muß die Größe D 
haben. Über die Aufteilung des Momentes D auf die einzelnen Zellen braucht jedoch nichts 
vorausgesetzt zu werden. Aus der Bedingung, daß die Arbeit des Zustandes mit den Schnitt- 
kräften N®, T® auf den Verschiebungswegen jeder einzelnen Konstanten CP verschwinden 
muß, ergeben sich bei einem Stabquerschnitt mit n Zellen » Gleichungen zur Bestimmung der 
n Unbekannten. Dabei zeigt sich, daß das entstehende Gleichungssystem bei den praktisch 
in Frage kommenden Querschnittsformen in dreigliedriger Form erhalten werden kann. 


Zur Ergänzung der bisherigen Untersuchungen über die torsionsfreie Verformung ist 
noch nachzuweisen, daß bei dem gewählten Spannungsansatz (5) die reine Zug- oder Druck- 
beanspruchung des Stabes durch eine Kraft P verdrehungsfrei ist. In (7) ist jetzt 


TU= TP=0 


zu setzen. Man erhält dann 


2 d Pteo 
BP 
2CcDP 
Psinwcosagdp. 


Mit sino= wird 


cos 


> ist, übereinstimmen und cos a 


da die obere und untere Grenze des Integrals, das gleich In (- 


stets positiv ist. 


8. Reine Verdrehung. Drillruhepunkt. 


Während die torsionsfreie Verformung des Stabes bereits zur Aufstellung des Ansatzes (5) 
klargestellt worden ist, ist von der reinen Verdrehbeanspruchung bisher nur der Spannungs- 
zustand, aber nicht der Verformungszustand betrachtet worden. Zur Ermittlung des letzteren 
sei wieder von der schon in Abschnitt 5 benutzten Annahme ausgegangen, daß durch die 
vorhandenen Versteifungen bzw. durch den Biegewiderstand der Mittelfläche der Schale eine 
Verformung vorgeschrieben wird, die der eines Stabes mit Vollquerschnitt entspricht. Die 
bekannte strenge Theorie der Torsion zylindrischer Stäbe liefert die Aussagen, daß die Ver- 
formung in einer Verdrehung der einzelnen Querschnitte besteht, wobei die Gestalt des Quer- 
schnittes in ihrer Projektion in Richtung der Stabachse exakt erhalten bleibt, während parallel 
zur Stabachse im allgemeinen eine Verwölbung auftritt. Es erscheint durchaus einleuchtend, 
das gleiche Verformungsbild auch bei stabförmigen Membranschalen vorauszusetzen, da z. B. 
aussteifende Rippen bei Beanspruchung senkrecht zu ihrer Ebene hinreichend weich sind, um 
eine Querschnittsverwölbung nicht merklich zu verhindern. Die in Frage stehende Annahme 
ist bei zylindrischen Membranschalen üblich und hat sich dort auch zweifellos bewährt ®). 
Überträgt man sie auch auf Kegelschalen, so wird die entstehende Lösung im allgemeinen 
nicht mehr wie bei Zylinderschalen mit den Verformungsbedingungen am einzelnen Schalen- 
element exakt verträglich sein®); die tatsächliche Verformung wird vielmehr nur im Mittel 


?) Über mehrzellige Querschnitte zylindrischer Stäbe vergleiche man W. J. Goodey: Airer. Engng. Bd. 8 
(1936), S. 93; A. Pflüger: a.a.0.; H. Ebner: Luftf.-Forschg. Bd. 14 (1937), S. 182; B. Radomski: Luftf.-Forschg. 
Bd. 17 (1940), S. 70. 

8) M. Vinsonneau: Rev. techn. Aeronaut. Bd 7 (1938), S. 7, empfiehlt allerdings, bei Berechnung der Ver- 
drehung zylindrischer Schalen gerade die entgegengesetzte Annahme behinderter Querscehnittsverwölbung an die Spitze 
der ganzen Untersuchung zu stellen. 

%) Ein Weg zur exakten Lösung des Torsionsproblems kegelförmiger Stäbe beliebiger Querschnittsgestalt wird 
gezeigt bei D. Panov: C. R. Acad. Sei. URSS., N. s. 20 (1938), S. 251. 
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eine Verdrehung darstellen. Die dadurch bei Kegeln neu hinzukommende Fehlermöglichkeit 
bewegt sich jedoch im Rahmen der Genauigkeit der übrigen Rechnung, da ja die gleiche Un- 
sicherheit für die Querkraft-Biegung auch schon bei zylindrischen Stäben besteht. 

Es sei y der Winkel, um den sich bei der Verformung ein Querschnitt aus seiner Null- 
lage heraus dreht. Die Ableitung von y nach « ist dann der gewöhnlich als Verdrehwinkel 
je Längeneinheit bezeichnete Winkel, um den sich zwei benachbarte Querschnitte gegen- 
einander verdrehen. Zur Durchführung der für die Ermittlung von y bzw. y’ erforderlichen 
Mittelbildung wird zweckmäßig die Formänderungsarbeit benutzt. Nach (7) ergibt sich die 
Arbeit, die der Zustand reiner Torsion mit dem Schubfluß 7? = =: beim Aufbringen 


des Drehmomentes D leistet, zu F 


1 G cos 
4 =3 Hits o) 


ode, 


wobei der Faktor 5 hinzugefügt wird, weil es sich um die Arbeit eines Spannungszustandes 


auf seinen eigenen Verschiebungswegen handelt. Ausgedrückt durch den Verdrehwinkel muß 
für die Formänderungsarbeit gelten 


1 
Durch Gleichsetzen der beiden Ausdrücke für A”? bekommt man 
D 1 „..,cos’o 
y $ 1+4tg’o) dep 


oder bezogen auf den Querschnitt an der Stelle «= a 
D a 


Nach der exakten Torsionstheorie zylindrischer Stäbe ist der Verdrehwinkel y’ eine für 


alle x konstante Größe. Man erkennt aus (10a), daß hier f= — sein muß, wenn (bei kon- 


stantem D) y" ebenfalls von « unabhängig werden soll. Wird die Verdrehsteifigkeit mit GJp 
bezeichnet, so folgt Jp aus der Definitionsgleichung G Jpy'’=D zu 


1 
$ 3 +4 tg’ ©) d 


Der Punkt, um den sich bei der Belastung durch ein Drehmoment D zwei benachbarte 
Querschnitte gegeneinander drehen, wird als „Drillruhepunkt“ bezeichnet. Da von der Wand- 
stärke ö vorausgesetzt war, daß sie sich in ihrem Verlauf längs der Querschnittskurve für 
verschiedene Querschnitte nur durch einen von « abhängigen Faktor unterscheidet, folgt, daß 
bei der Kegelschale die Verbindungslinie der Drillruhepunkte, die Drillachse, eine durch die 
Kegelspitze hindurchgehende Gerade sein muß. Im folgenden sei 


— Abstand des Drillruhepunktes von der 2-Achse, 


u Verschiebung eines Punktes der Schalenmittelfläche in Richtung der Erzeugenden 
bei reiner Verdrehbeanspruchung. 


Zur Ermittlung des Drillruhepunktes werden die Ergebnisse des Abschnittes 7 über torsions- 
freie Verformung benutzt. Zunächst sei kurz auf den einfachen Fall eines zylindrischen 
Stabes eingegangen. Hierbei ergibt sich bekanntlich, daß die Querkraftmittellinie und die 
Drillachse identisch sind. Zum Beweise wird gewöhnlich ein eingespannter Balken (Bild 1) 
betrachtet‘), der an seinem freien Ende nur durch eine in Richtung der z-Achse wirkende 
Einzellast beansprucht wird. Der Angriffspunkt der Einzellast soll der Querkraftmittelpunkt 
sein, so daß die Arbeit zu Null wird, die der im Stab durch die Belastung hervorgerufene 
Spannungszustand bei einer reinen Verdrehung leistet. Diese Bedingung läßt sich aber nicht 


10) Vgl. z. B. A. Föppl: Drang und Zwang Bd. 2, München und Berlin 1928, S. 129. 
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nur durch die Arbeit der inneren Kräfte, sondern auch durch die der äußeren Kräfte aus- 
drücken. Die äußeren Kräfte bestehen am freien Ende des Stabes aus der Einzelkraft und 
an der Einspannstelle aus der Auflagerkraft und dem Einspannmoment. Die Auflagerkraft 
leistet bei einer Verdrehung des Stabes von vornherein keine Arbeit. Die Arbeit der Einzel- 
last wird dann zu Null, wenn sich der Stab gerade um die Querkraftmittellinie dreht. Um 
aussagen zu können, daß Drillachse und Querkraftmittellinie zusammenfallen, muß also nur 
noch bestätigt werden, daß die Arbeit verschwindet, die das Einspannmoment auf den Ver- 
schiebungswegen der zur Torsion gehörigen Verwölbung leistet. Die Erfüllung dieser Bedingung 
läßt sich in der Theorie der Torsion zylindrischer Stäbe bei Berechnung der Verwölbungen u 
leicht nachweisen ''). Mit den hier benutzten Bezeichnungen lautet die Bedingung für eine 
zylindrische Membranschale 


Pöuzody=V 


Die Betrachtung sei jetzt auch auf Kegelschalen ausgedehnt. Hier wird man nicht mehr 
voraussetzen dürfen, daß eine Gl. (11) entsprechende Beziehung uneingeschränkt Gültigkeit 
hat. Es läßt sich nämlich schon unmittelbar durch die Anschauung bestätigen, daß die Punkte 
einer Querschnittskurve auch dann Verschiebungen « erleiden, wenn sich ein Stabelement, 
wie es etwa in Abb. 8 dargestellt ist, als starrer Körper um die Drillachse dreht, d. h. also, 
wenn y’—=0 aber y=F0 ist. Bei einer derartigen Verwölbung, die eine Drehung der Quer- 
schnittsebene um die y-Achse mit sich bringt, wird aber zweifellos ein in dem betreffenden 
Querschnitt angreifendes Biegemoment Arbeit leisten. Es 
dürfte jedoch durchaus plausibel sein, wenn man in Über- a 
einstimmung mit den bisherigen Annahmen über die Ähn- ! y 
lichkeit der Verformungen eines kegelförmigen Stabes und 
eines gewöhnlichen zylindrischen Stabes voraussetzt, daß 
für y=0, also z. B. an der Einspannstelle eines Frei- 
trägers die Verschiebungen « derart beschaffen sind, daß 
mit hinreichender Näherung die Arbeit eines Biege- ) / 


Bild 8. Stabelement. 


Stabochse 


momentes bei der Verwölbung zu Null wird. Es soll ll Brstsches 
also angenommen werden F An 


Die Lage des Drillruhepunktes ergibt sich nun, wenn man die 
bei Verdrehung geleistete Arbeit eines Stabelementes gleich 
Null setzt, das nach Bild 8 durch Biegemomente und durch 

Querkraftmittelpunkt angreifende Querkräfte belastet 
wird. Die durch d. dividierte Arbeit des Stabelementes ist = erde 


cos a 


Das letzte der beiden Integrale kann, wie folgt, umgeformt werden: 
| 
j 2.) udpo= \cos [M.+9 (x udg 


cos a Say? 


- zuody. 
cos a 


M 


-zwWod ——zvodpo=0. 
cos a +7 cos a 


M 
= elle 


Man erhält dann 
Q [fe — ec) —e' vj+7 


Die so gewonnene Beziehung muß für jeden FE Wert des Drehwinkels y und 
des Biegemomentes M gelten. Insbesondere muß sie auch dann richtig sein, wenn infolge 
geeigneter Stützung und Belastung des Stabes an der betrachteten Stelle x die Größen y 
und M gerade zu Null werden. In diesem Fall gilt dann 


11) Vgl. z. B. H. Neuber: Z. angew. Math, Mech. Bd. 21 (1941), S. 92, 
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0) 
(e— ey=0) y +7 


Daraus folgt bei Benutzung der Annalıme (12) 


und nach (9) und (9a) 


1 
0 


Die Koordinaten des Drillruhepunktes und des Querkraftmittelpunktes stimmen also dann 
überein, wenn der Stab an der betreffenden Stelle nur durch eine Querkraft, aber nicht 
durch ein Biegemoment beansprucht wird. Die Ausdrücke (9) und (9a) können damit unter 
Benutzung von (13) und (13a) in der kürzeren Form 


(Ca — 23 Sa a) + 2 Fa 
geschrieben werden. 

Ist @ sehr klein, so wird die Benutzung der Größe ce für praktische Rechnungen 
unzweckmäßig. Zur Ermittlung von D kann dann folgende Beziehung von Interesse sein. 
Es sei mit D, die Summe der Momente aller in einem Querschnitt übertragenen Schnittkräfte 
in bezug auf den Drillruhepunkt bezeichnet. Zwischen D, und der bisher benutzten Größe D 
besteht der Zusammenhang 


D=D,+9(e—e) 
=D. +9e— Yle—2%k, x) — 


M 
=D, 


Für mehrzellige Querschnitte behalten alle für die reine Verdrehbeanspruchung bereits 
gewonnenen Ergebnisse weitgehend ihre Gültigkeit. Die durchzuführende Rechnung erschwert 
sich lediglich dadurch, daß das Drehmoment D auf die einzelnen Zellen den getroffenen An- 
nahmen entsprechend aufgeteilt werden muß. Da vorausgesetzt war, daß die Stabverformung 
bei Torsionsbeanspruchung in einer Verdrehung besteht, muß die Aufteilung so erfolgen, daß 
sich alle Zellen um den gleichen Winkel y’ verdrehen. 


9. Vereinfachungen für die praktische Rechnung. 

In den vorstehenden Abschnitten ist eine ganze Reihe von Annahmen enthalten, die der 
technischen Balkenbiegungslehre entnommen sind. Die besonderen durch die stabförmige 
Gestalt des Kegels bedingten Beziehungen für die Spannungsverteilung am einzelnen Schalen- 
element sind dagegen ohne irgendwelche Vernachlässigungen aufgestellt worden. So notwendig 
diese Exaktheit zunächst war, so ergibt sich doch nachträglich die Möglichkeit, die abgeleiteten 
Formeln für praktische Rechnungen noch zu vereinfachen. Infolge der stabförmigen Gestalt 
des untersuchten Kegels werden die Winkel a und o, die für einen Zylinder zu Null werden, 
immer kleine Größen sein. Es wird daher häufig zulässig sein, cos « und cos & angenähert 
gleich Eins zu setzen und sin a bzw. sin gegenüber Eins zu vernachlässigen. Es wird dann 


F-F, 
.Damit wird aus den Gl. (5) und (6) 


r 


Der Ansatz für die Normalspannungen ist jetzt der gleiche wie bei einem zylindrischen Stab; 
bei den Schubkräften sind jedoch die besonderen Vorschriften über die Berechnung von Ü 
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und die Abminderung der Querkraft nach wie vor zu beachten. Für die Größen k, und k, 
bekommt man 


ode 


Für die Lage des Drillruhepunktes ergibt sich damit die Formel 


Pl özedp)redg, 


die wieder mit der für einen Zylinder gültigen Beziehung übereinstimmt. Für den Querkraft- 
mittelpunkt bekommt man 


Da ec nicht mehr von x abhängig ist, muß die Querkraftmittellinie eine Parallele zur Stabachse 
sein. Das bedeutet, daß man jetzt für den ganzen Stab und nicht nur für eine bestimmte Stelle 
bei geeigneter Belastung torsionsfreie Verformung erreichen kann. Die Beziehung für ce dürfte 
anschaulich besonders einleuchtend sein: Da das Moment von 9 um die Stabachse gleich 
dem Produkt Ge ist, muß der Angriffspunkt des Quer- 

Drillnuhepunkt Querkroftmittelgunkt  kraftanteils der Drillruhepunkt sein. Man bekommt 


r also die in Bild 9 dargestellten Verhältnisse. Die 
Resultierende der Schubkräfte wirkt in demselben 
4 / Punkt, durch den sie auch bei einem zylindrischen 


2 Stab hindurchgeht. Der ganze Unterschied in der 
_. Berechnung eines Kegel- und eines Zylinderquerschnitts 


Bild 9. Angriffspunkte der Querkraftanteile bei jst nur noch durch den im Schwerpunkt angreifenden 
torsionsfreier Verformung nach der Näherungs- 


rechnung. Querkraftanteil = bedingt. Zwischen dem Dreh- 


moment D um den Querkraftmittelpunkt und dem Drehmoment D, um den Drillruhepunkt 
besteht der Zusammenhang 


M 
DA D.+7e. 


Unter Berücksichtigung dieser Beziehung läßt sich nun die gesamte durchzuführende Rechnung 
in der folgenden einfachen Regel zusammenfassen: 


Die Spannungsberechnung einer stabförmigen Membran-Kegelschale kann in erster 
Näherung genau so wie für einen zylindrischen Stab erfolgen; es ist lediglich an 


M 
Stelle der Querkraft 9 eine um den Betrag > verminderte Querkraft und an 


Stelle des Drehmomentes um den Drillruhepunkt ein um den Betrag = e ver- 
größertes Drehmoment einzusetzen. 


10. Versuche mit einer Kegelschale. 


Zur Nachprüfung der in den vorstehenden Abschnitten entwickelten Theorie wurden an 
einer Kegelschale aus Leichtmetall (Al-Cu-Mg-Legierung) Spannungs- und Verformungsmessungen 
vorgenommen '*). Die Abmessungen des Versuchsstückes gehen aus Bild 10 hervor. An beiden 
Enden ist die Schale durch kräftige Stahlplatten abgeschlossen, während sie im Innern durch 
gebördelte Rippen ausgesteift wird. Die Querschnittskurve hat rechteckige Gestalt. Die 
Wandstärke beträgt für drei Seiten des Rechtecks 3 mm und für die vierte (kurze) Seite 1 mm. 
Es läßt sich nachweisen '°), daß eine zylindrische Schale gleicher Querschnittsabmessungen 
verwölbungsfrei ist. Nach den Ausführungen von Abschnitt 8 wird dann bei der Kegelschale 
die Verwölbung mit guter Näherung so vor sich gehen, daß die Querschnitisebene bei der 
Verdrehung eine Ebene bleibt. Die die Schale abschließenden Stahlplatten werden also durch 
Verwölbungskräfte nicht beansprucht, so daß eine Differenz zwischen Theorie und Versuch 
infolge behinderter Querschnittsverwölbung nicht befürchtet zu werden braucht. Die in den 


12) Diese Versuche wurden in der Werkstoffversuchsabteilung der Focke-Wulf Flugzeugbau G. m.b. H., Bremen, 
von H. Mester durchgeführt. 
13) H. Neuber: Z. angew. Math. Mech. Bd. 21 (1941), S. 9. 
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Ecken der Schale zur Verbindung der Bleche erforderlichen Winkel und die zur Verstärkung 
der Einspannstelle aufgenieteten Laschen bedingen allerdings eine geringe Änderung der 
Verhältnisse, die jedoch, wie die Versuchsergebnisse zeigen, hinreichend klein ist, um sich 
für die Spannungsverteilung im mittleren Teil der Schale nieht mehr bemerkbar zu machen. 
Der Kegel wurde an dem Endquerschnitt mit den kleineren Abmessungen an einem Gerüst 

fest eingespannt, während er an dem 

anderen Ende belastet wurde. Die 
L 30x30x3 gewählte Art der Lagerung und Kraft- 
einleitung ist aus dem Lichtbild, 
Bild 11, zu erkennen, das die Schale 
bei Verdrehbeanspruchung zeigt. Die 
u a Einleitung der Kräfte erfolgte am 
Endquerschnitt mit den größeren Ab- 
messungen, weil auf diese Weise der 
Angriffspunkt einer Kraft im Ver- 
hältnis zu den Querschnittsabmessun- 
gen am genauesten festgelegt werden 
konnte. 


* w Stellen, an denen Dehnungsmessungen durchgeführt wurden 
Bild 10 (oben). Skizze des Versuchsstückes. 


Bild 11 (rechts). Lichtbild des Versuchs- 
stückes bei Verdrehbeanspruchung. 


Zur Ermittlung des Spannungs- 
zustandes wurden für einen (uer- 
schnitt zwischen den beiden Rippen 
Dehnungsmessungen an den in Bild 10 
angegebenen Stellen durchgeführt. 
(Vgl. auch Bild 11, auf der einige 
Meßgeräte zu erkennen sind.) Obwohl 
der Querschnitt eine Symmetrieachse 
besitzt, wurden doch in beiden Quer- 
schnittshälften Meßstellen angeordnet, 
um die Möglichkeit zu haben, durch 
Mittelbildung der sich entsprechenden Meßwerte den Einfluß von Ungenauigkeiten in der 
Form der Schale ausgleichen zu können. Die in der Theorie enthaltene Voraussetzung eines 
Membranspannungszustandes verlangt, daß Biegespannungen in der Schalenfläche nicht 
auftreten können. Praktisch sind jedoch stets derartige Spannungen vorhanden. Für die 
Dimensionierung der Schale sind sie zwar ohne Interesse, da das Material bei Überschreitung 
der Fließgrenze in den äußeren Fasern des Bleches durch Spannungsausgleich von selbst 
dafür sorgt, daß die von der Membrantheorie gelieferten Mittelwerte der Spannungen ein- 
gehalten werden; bei Spannungsmessungen im elastischen Bereich müssen sie jedoch be- 
rücksichtigt werden, wenn man eine brauchbare Übereinstimmung zwischen Theorie und 
Messung erhalten will. Das übliche Verfahren zur Elimination der Biegespannungen besteht 
darin, daß an einer Stelle der Schale der Spannungszustand auf beiden Seiten des Bleches 
gemessen und von den Messungen das Mittel gebildet wird. Diese Methode wär im vorliegenden 
Fall nicht anwendbar, da das Versuchsstück nur von einer Seite zugänglich war. Die Er- 
mittlung der Biegungsspannungen war jedoch durch Verwendung besonderer Meßgeräte möglich, 
die außer der Dehnung auch noch die Beulung des Bleches anzeigen '*). 


In den Bildern 12 und 13 sind bei zwei verschiedenen Beanspruchungszuständen für 
die Längskräfte N und die Schubkräfte 7 die Ergebnisse von Theorie und Messung einander 
gegenübergestellt. Die beiden Spannungszustände wurden dadurch erzeugt, daß am freien 
Ende der Schale eine Einzelkraft von 1000 kg einmal in Richtung der y- und einmal in 
Richtung der 2-Achse wirkte. Die Übereinstimmung zwischen Theorie und Versuch kann als 
durchaus befriedigend bezeichnet werden. Die in den Bildern 12b und 13b eingetragenen 
gestrichelten Kurven geben den Verlauf der Schubkräfte 7 wieder, wie er sich nach einer 
Rechnung ergibt, bei der der Kegel als zylindrischer Stab mit den Abmessungen des Me&- 


14) Es handelt sich hierbei um ein Meßinstrument, das in der Werkstoffversuchsabteilung der Focke-Wulf 
Flugzeugbau G. m. b. H., Bremen, entwickelt wurde. Näheres s. h. A. Dose, Luftf.-Forschg. Bd. 18 (1941), S. 95. 
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querschnitts angesehen wird. Man erkennt, daß der Fehler, den man auf diese Weise begeht, 
recht beträchtlich ist. Der Unterschied zwischen den Ergebnissen der in den Abschnitten 
1 bis 8 entwickelten genaueren Rechnung und der Näherungsrechnung des Abschnittes 9 ist 
für die in den Bildern 12 und 13 angegebenen Spannungszustände so gering, daß er sich in 
dem gewählten Maßstab nicht darstellen läßt. Für die Lage des Querkraftmittelpunktes 
ergibt sich jedoch ein etwas größerer Unterschied. Nach der Näherungsrechnung würde für 
den Meßquerschnitt 


— 39 [mm] 


sein, während nach der genaueren Rechnung 


=— 17-5 +8 [mm] 


Bild 12a (links). Längskräfte 

N bei Beanspruchung aus 

Querkraft und Biegung in 
der x-y-Ebene. 


Bild 12b (rechts). Schub- 
kräfte T bei Beanspruchung 0 #0 Rgjem 
Kröftemaßstob: 0 0 60 Bokglım aus Querkraft und Aröfteverlauf m Rechn 
wii J U 
Kröftevenlauf Rechnung “in der x-y-Ebene. © Meßounkte 


Meßpunkte für eme zylindr Schole 


610mkg 


Bild 13a (links). Längskräfte 

N bei Beanspruchung aus 

Querkraft und Biegung pa- 
rallel zur x-z-Ebene. 


700049 
= 


Bild 13b (rechts). Schub- 0 20 40 60 80kg/cm 
kräfte T bei Beanspruchung #röftemaßstab:_ 11 


Kräftererlauf n. Rechnung aus Querkraft \und „Biegung Aröfteverlauf n Rechnung 
© Meßounkte parallel zur x-z-Ebene. e Meßpunkte 
— — Aröfteverlauf für eıne zylinde Schale 


; ist. In Bild 14 sind Messungen über den bei reiner 
RR Torsion auftretenden Drehwinkel y den gerechneten 
_ Werten gegenübergestellt. Die Verformung wurde dadurch 
Drillachse Ai j gemessen, daß an der Schale Maßstäbe befestigt wurden, 
en deren Verschiebung mit Hilfe eines Nivellierinstrumentes 

beobachtet wurde. Nach der Theorie erhält man für y 
die Gleichung 


1,05 
y’=133 — (x in m). 


Sie gilt nicht mehr in der Nähe der Einspannstelle im Be- 
reich der dort vorhandenen Verstärkungslaschen. Dieser 
® Mebpunkte Bereich wurde zu einem Vergleich zwischen Rechnung und 
m Messung nicht herangezogen. Die Integrationskonstante 
r in der Gleichung für y wurde so gewählt, daß an der ' 

Bild 14. Drehwinkel bei einem Drehmoment Y 
von 500 mkg. Rippe bei «= 1,2m Versuch und Rechnung übereinstimmen. 


11. Zusammenfassung. 


Bei der Spannungsberechnung von Stäben, die einen dünnwandigen Hohlquerschnitt 
mit längs der Stabachse veränderlichen Abmessungen besitzen, ist es im allgemeinen nicht 
mehr ausreichend genau, wenn man den zu untersuchenden Stab an der gerade betrachteten 
Stelle durch einen Zylinder annähert. Es wird vielmehr erforderlich, in der Näherung einen 
Schritt weiter zu gehen, und den Stab durch eine Kegelschale zu ersetzen. Die wichtigsten 
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Ergebnisse, die in den Abschnitten 1 bis 9 bei Übertragung der Ansätze der technischen 
Balkenbiegungslehre auf kegelförmige Schalen erhalten wurden und, wie Abschnitt 10 zeigt, 
durch Dehnungsmessungen an einem Versuchsstück bestätigt werden konnten, sind die folgenden. 
Die in der Schale in Richtung der Erzeugenden auftretenden Längskräfte können — von 
kleinen praktisch meist nicht wesentlichen Korrekturgliedern abgesehen — wie bei einem 
Zylinder berechnet werden. Für die Schubspannungen ergeben sich jedoch erhebliche 
Änderungen, die in der Hauptsache darin bestehen, daß statt der Querkraft und statt des 
Drehmomentes eines zylindrischen Stabes andere Größen maßgeblich werden, die sich von 
den ersteren durch einen vom Biegemoment abhängigen Betrag unterscheiden. Der Drill- 
ruhepunkt eines Kegelquerschnittes stimmt im wesentlichen mit dem Drillruhepunkt eines 
entsprechenden Zylinderquerschnitts überein. Die Lage des Querkraftmittelpunktes ist jedoch 
von der Belastung abhängig. Drillruhepunkt und Querkraftmittelpunkt fallen lediglich dann 
zusammen, wenn der betreffende Stabquerschnitt nur durch eine Querkraft, aber nicht durch 
ein Biegemoment beansprucht wird. Bei Belastung allein durch ein Biegemoment ist die 
Verformung des Stabes im allgemeinen nicht torsionsfrei. 305 


KLEINE MITTEILUNGEN 


Fehlerbetrachtung zur Rinaldinischen 
Regel. Nach der Rinaldinischen Regel kann ein 
Kreis in beliebig viele näherungsweise gleiche 
Teile geteilt werden, und im folgenden wird die 
Größe der hierbei auftretenden Teilwinkelfehler 
näher untersucht. 


Zur Konstruktion eines Sehnen-n-eckes wird der 


hierin ist 


(aus ACMP), 


ferner 
2-/r-cosa (aus ACDM) 


und 


Durchmesser in n gleiche Teile geteilt und von 
den Schnittpunkten O und P (der Kreisbögen um 
die Durchmesserendpunkte A und B mit AB als 
Halbmesser) die Strahlen über jeden zweiten Durch- 
messerteilpunkt gezogen, die dann den Umfang des 
Kreises in den gesuchten Teilpunkten treffen 


(Bild 1). 
A b 


Setzt man nun für 


—=Yd?+e?— de 


Durch Einsetzen und Umstellen erhält man 


worin m die Ordnungszahl der Teilung und n die 
Seitenzahl des Vieleckes ist, so ergibt sich nach 
Kürzen von r 


Nun führt man die Hilfsgröße 


n 


ein und erhält 


Bild 2. 


z =| 3+ 2? 


Zur Ermittlung des Teilwinkelfehlers kann der Teil- 
winkel aus den Konstruktionsgrößen bestimmt 


werden und zwar ist (Bild 2) eine quadratische Gleichung mit der Lösung 


3 + 922 


3-+ 2? 


(b+ in sind 


cosa—sin (MV '+a) = (aus A DMP), 


(aus A ACP). 


Winkelfehler 


— 
- 


Winkelfehler 
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wobei das negative Wurzelvorzeichen für den 
Schnittpunkt E des Strahles mit dem Kreisumfang 
gilt, und daher für die Fehlerrechnung nicht be- 
rücksichtigt zu m.. braucht. Da der Winkel bei 
fehlerfreier Teilung - 360 ° betragen soll, ergibt 


sich schließlich 


cos (da+" 360°) = 


Wegen der Symmetrie der Konstruktion ist die 
Fehlerermittlung nur für 0°<Sa< erforderlich, 
d. h. also für 
in Bild 3 zusammengestellt. 


Die Fehlerwerte sind 


Teilwinkel 
60° 65° 70° 75° 80° 85° 90° 
5° m° 15° 20° 25° 30° 35° 40° 45° 50° 55° 


Bild 3. Fehlerschaubild zur Rinaldinischen Regel. 


Man ersieht hieraus, daß der Fehler etwa bei 18° 
am größten ist und rund 38° beträgt. Bemerkens- 
wert ist der geringe Fehler im Gebiet zwischen 
60° und 90°. Für Konstruktionen, bei denen ein 
Winkelfehler bis zu 3° zulässig ist, was in der 
Praxis sehr häufig der Fall sein wird, kann also 
die Rinaldinische Regel durchaus angewendet wer- 
den, wenn man zur Ermittlung der Sehnenlänge 
nur Teilpunkte im Bereich der kleinen Fehler wählt. 


Hamburg. H. Peineke. 328 


über Strömungen von Schäumen. Unter 
einem Schaum soll im folgenden das innige Ge- 
misch einer Flüssigkeit und eines Gases verstanden 
werden. Es soll keine merkliche Entmischungs- 
tendenz vorhanden sein und kein nennenswerter 
Übergang der flüssigen Komponente in die gas- 
förmige oder umgekehrt stattfinden. Die Mischung 
soll so innig sein, daß Volumselemente, deren Ab- 
messungen genüige end klein sind gegenüber den Ab- 
messungen der Strömung, noch immer Mischungs- 
charakter besitzen. 

Praktische Bedeutung besitzen solche Strömungen 
vor allem bei Wasserstrahlluftpumpen. Theoretisch 
sind sie noch wenig behandelt. Ansätze finden 
sich bei Pfleiderer t). 

Es sollen zunächst die hydrodynamischen Grund- 
gleichungen der Schaumströmung entwickelt wer- 
den, unter der Voraussetzung, daß die flüssige und 
gasförmige Komponente in jedem Volumselement 
in demselben Gewichtsverhältnis steht. Greift man 
ein Volumen V des Schaumes heraus, so sollen 
V; und V, sowie G; und Gy die Volumsanteile 
sowie die Gewichtsanteile der flüssigen und gas- 
förmigen Komponente bedeuten. Es ist dann das 
Mischungsverhältnis: 


G 


Für eine vorliegende Strömung nach Voraussetzung 
von V unabhängig und konstant. Es sei ferner die 
gasförmige Komponente ein ideales Gas, die flüs- 


Zur Berechnung der 
luftpumpe. Z. VDI Bd. 58 (1914) S. 956 u. S. 5 


sige Komponente sei unzusammendrückbar. Unter 
diesen Voraussetzungen ist die Zustandsgleichung 
des Schaumes leicht abzuleiten. Bezieht man die 
eingeführten Größen Gg, Gy, Vg. Vz; auf das Volums- 
element und sind og, o und e die Dichten der gas- 
förmigen und der flüssigen Komponente sowie des 
Schaumes, so gilt: 


Gr+Gs 


1’ 


ferner gilt, wenn ? der Druck, R die Gaskon- 
stante und 7 die absolute Temperatur bedeutet: 


p 
(2) 


Daraus folgt: 
1 1 Br 

Pe 
als Zustandsgleichung des Schaumes, worin or als 
konstant angesehen werden kann, wenn im 
Temperaturbereich, der in Frage kommt, die 
Dichteänderungen der flüssigen Komponente be- 
langlos sind. 

Setzt man vollkommene Konvektion voraus, so 
daß die flüssige und gasförmige Komponente an 
jeder Stelle der Strömung eine einheitliche Ge- 
schwindigkeit besitzen, so gelten, wenn man von 
Reibungsverlusten absieht und wenn die äußeren 
Kräfte ein Potential U besitzen, die Eulerschen 
Gleichungen und die Kontinuitätsgleichung wie 
bei der Strömung des Kontinuums. Es gilt also, 
in vektorieller Form geschrieben ?): 


. (4), 


d0 
. (4a), 


wobei % den Geschwindigkeitsvektor der Schaum- 
strömung bedeutet. 


Um den ersten Hauptsatz dq—=du+ Apdv ein- 
zuführen, muß noch der Ausdruck für die innere 
Energie du gefunden werden. Es ist du, bezogen 
auf die Gewichtseinheit des Schaumes: 


du=|c, 
| ) 


Gy 
wenn €, die spezifische Wärme der gasförmigen 
Komponente bei konstantem Volumen und cr die 
spezifische Wärme der flüssigen Komponente be- 
zeichnen. Nach (1) ist also: 


und damit der erste Hauptsatz: 


1 


wobei die Beziehung o-., zu beachten ist. 


Um ein geschlossenes Gleichungssystem zu er- 
halten, fehlt noch eine Gleichung, die den Wärme- 
transport festlegt. Es wird für praktische Zwecke 
hinreichen, wenn man den Wärmeübergang ganz 

2) Es bedeutet @5 das Innenprodukt, xD das Vektor- 
produkt und @,5 das dyadische Produkt. 
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vernachlässigt. Man erhält dann eine adiabatische 
Zustandsänderung, die durch dq=0 bestimmt ist. 
Führt man dies in (5) ein, so erhält man, bei Be- 
nutzung von (3), die Differentialgleichung der 
Adiabate: 


dp do 
P 
or (1-+ u) 


wobei von der bekannten Beziehung (»—=A'R—c,. 
Gebrauch gemacht wurde. Die Integration von (6) 
liefert die Adiabatengleichung des Schaumes: 


(6), 


1 1 . . 


oe 
Für die beiden Grenzfälle «—0 (reine Flüssigkeit) 
und «=» (reines Gas) erhält man hieraus die 


bekannten Beziehungen = eonst und gr 
c . . 
eonst, wenn 2#= > bedeutet. Die Gleichungen (4), 


(4a) und (7) reichen zur eindeutigen Beschreibung 
der Strömung hin, falls die nötigen Randwerte 
festliegen. 


Wesentlich für die Beurteilung des Strömungs- 
charakters ist die Frage, ob sich die Strömung im 
unterkritischen oder überkritischen Bereich be- 
findet. Die Grenze zwischen beiden Bereichen liegt 
bekanntlich dort, wo die Strömungsgeschwindigkeit 
w mit der Schallgeschwindigkeit übereinstimmt. 
Diese ist nach einer bekannten Beziehung gegeben 


Berechnet man 
do do 


ergibt sich in Verbindung mit (2): 


durch @ = aus (6), so 


1 BTV_ 
(men te)-gRT 
woraus man entnehmen kann, daß die Schall- 


geschwindigkeit des Schaumes im Gegensatz zu 
der des idealen Gases druckabhängig ist. 


(8), 


Führt man als dimensionslose Variable: 


o-gRT' 


ein, so erhält man aus (8): 


(1-+u) (u+») 


a 
VrgRT 


wobei: 


=(P+ 


(N) 
Cy 


bedeutet. 


Für den Grenzfall der reinen Flüssigkeit (« = 0) 
wird @—=x, da die Flüssigkeit als unzusammen- 
drückbar vorausgesetzt wurde. Für den Grenzfall 
des reinen Gases («= &) erhält man die bekannte 
Gleichung e=yx9gRT. Aus (10) ergibt sich aber 
die bemerkenswerte Tatsache, daß für gewisse Be- 


reiche des Mischungsverhältnisses « die Schall- 
geschwindigkeit des Schaumes wesentlich Kleiner 
werden kann als die der gasförmigen Komponente 
bei gleicher Temperatur. 


Faßt man nach (10) a@ bei konstantem P und T 
als Funktion von « auf, so besitzt a für einen be- 
stimmten Wert wein Minimum. Für ı folgt aus (10) 
die Beziehung: 


Aus (9) und (2) folgt ferner: 


p—® 


(9a), 


woraus man ersieht, daß man P< 1 voraussetzen 
darf. Da x und v von der Größenordnung eins 
sind, „kann man statt (12) angenähert schreiben: 


2xv v? 

3 

= v? 

—v 


(12a). 


Um von der Größenordnung der Koeffizienten ein 
Bild zu bekommen, sollen die Werte für Wasser 
und Luft eingesetzt werden #=5'8, 14): 

+ + 876R 


Da P<1 ist, existiert nur eine positive Wurzel, 
die sich näherungsweise zu: 
(13) 


ergibt. Daraus folgt für Ayin: 


. . (10a) 
P(1+P)(P+») 


oder näherungsweise: 


VzgRT % 


Man sieht aus (10b), daß “in nur einen kleinen 
Bruchteil der Schallgeschwindigkeit der gasförmigen 
Komponente beträgt. 


(10b). 


Um die Bedingung (13) anschaulich zu deuten, 
ersetzt man sie nach (1) und (9a) durch: 


Ga _ 00 

Gr 
daraus folgt: Die Schallgeschwindigkeit 
des Schaumes ist also dort ein Minimum, wo die 
flüssige und gasförmige Komponente je das halbe 
Volumen eines Elementes einnehmen. 
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Buchbesprechungen 


BUCHBESPRECHUNGEN 


L. CESARI, F. CONFORTO und C. MINELLI, 
Travi Continue Inflesse e Sollecitate 
Assialmente. Mit einem Vorwort von 
A. Fiore. (Publicazioni Dell’Istituto Per Le 
Applicazioni Del Calcolo, Consiglio Nazionale Delle 
Ricerche Roma 1941/XIX, N. 91.) VII+244 S. 
Rom 1941. Preis geb. 65 Lire. 

Die Dreimomentengleichungen von G. Albenga 
— die bekannten Verallgemeinerungen der Ula- 
peyronschen — für durchlaufende Balken, die in 
der Querrichtung auf Biegung und gleichzeitig in 
der Längsrichtung auf Zug oder Druck (P) bean- 
sprucht werden, bilden die Grundlage für die Be- 
rechnung der Flugzeugholme für diagonal ver- 
spannte Doppeldecker; dabei genügt es in den 
meisten Fällen, die Trägheitsmomente (J) der 
Holmquerschnitte für die einzelnen Felder kon- 
stant anzunehmen. Schon bei Eindeckern führt 
jedoch die elementare Bedingung gleicher Biege- 
steifheit auf eine Veränderlichkeit der Holmhöhe 
und damit von J. Obwohl diese Konstruktionen 
ihre vormalige Bedeutung in weitem Maße ver- 
loren haben, werden diese Gleichungen immer ihre 
Bedeutung behalten. Die Verfasser geben in der 
vorliegenden Darstellung eine übersichtliche und 
durch Zahlen- und Kurventafeln ergänzte Dar- 
stellung des mit diesen Gleichungen verbundenen 
Aufgabenkreises. Die in den Dreimomentengleichun- 
gen von G. Albenga auftretenden Konstanten 
hängen unter den angegebenen Bedingungen von 
zwei Parametern ab, nämlich von «= (Jyar — 
Imind/Imin und «= Pl2/EJ „in, so daß die Zahlen- 
tafeln einen doppelten Eingang erhalten. — Zum 
Schlusse wird ausführlich die vollständige analy- 
tische Lösung des Problems mit linear veränder- 
lichem Trägheitsmoment gegeben, die durch Bes- 
selsche Funktionen geleistet wird. 
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Dr. H. HAPPEL, o. Prof. a. d. Technischen Hoch- 
schule Breslau, Das Dreikörperproblem. 
Vorlesungen über Himmelsmechanik. 
XI + 526 S. m. 17 Fig. Leipzig 1941. K. F. Koehler 
Verlag. Preis geb. 30 M. e 

Es ist auch vom Standpunkt des angewandten 
Mathematikers aus sehr zu begrüßen, daß ein Buch 
erscheint, das den Studierenden in die Probleme 
der Himmelsmechanik und insbesondere in das 
Dreikörperproblem einführen will. Bekommt doch 
dadurch die allgemeine Mechanik den altklassi- 
schen Abschluß und kann auch der Nichtastronom 
daraus für andere Zwecke die Integrationsmethoden 
kennen lernen, die die Astronomen im Laufe einer 
langen Entwicklung ausgebildet haben. So kann 
von diesem Buche eine allgemeine, über das be- 
sondere Fachgebiet hinausgehende Befruchtung der 
Integrationstheorie der Bewegungsgleichungen der 
Mechanik ausgehen. Aus diesem Grunde ist es 
sehr erfreulich, daß das Buch die Grundlagen der 
klassischen Mechanik, Lagrangesche Gleichungen, 
kanonische Gleichungen und ihre Transformations- 
theorie in den beiden ersten Kapiteln von Anfang 
an entwickelt. Das dritte Kapitel beginnt dann 
das Hauptthema mit dem Zweikörperproblem, im 
vierten Kapitel werden die Delaunayschen und 
Poincar6schen Koordinaten für die gestörten 
Bahnelemente eingeführt, im fünften Kapitel die 
Störungen der Planeten behandelt. Das sechste 
und siebente Kapitel über den Mond enthält eigene 
Untersuchungen des Verf. Im Kapitel VIII folgen 
die periodischen Lösungen nach Methoden von 
Poincare, Wilkens und E. Hopf, in Kapitel IX 
Poincar6es Theorie der charakteristischen Ex- 
ponenten. Kapitel X behandelt weitere Fälle, in 


denen eine strenge oder angenäherte Interration 
möglich ist, auch die Frage der Kommensurabili- 
täten. Das letzte Kapitel geht noch kurz auf die 
Regularisierung nach Thiele, Levi-Civitä und 
Sundmann ein. Benutzt wird dabei die Birkhoff- 
sche Darstellung. Ein abschließendes Urteil über 
dieses sicher reichhaltige und anregende Buch 
muß dem Astronomen von Fach vorbehalten bleiben. 
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LEIF JOHNSEN, Dynamique Göncrale 
des systemes non-hologomes. Schriften, 
herausgegeben von der norwegischen Akademie 
der Wissenschaften in Oslo, mathematisch-natur- 
wissenschaftliche Klasse, 1941, Nr. 4. Oslo 1941, 
Verlag in Kommission Hos Jacob Dybwad. 75 
S. Preis 6,— Kr. 

Es ist zu begrüßen, daß der Verf. seine wert- 
vollen Untersuchungen über nichtholonome Systeme 
in umfassenderer Form herausgegeben hat. Im 
ersten Kapitel bespricht er neben der allgemeinen 
Einführung die Methoden des Ref., ferner diejeni- 
gen von Woronetz, Appell und Schieldrop. Es 
werden gleich von vornherein auch nichtlineare 
Bedingungen zugelassen. Wenn Ref. hier gleich 
eine kritische Bemerkung machen darf, so die, daß 
die Methode von ihm und die von Woronetz beide 
in einer ganz allgemeinen enthalten sind, die Ref. 
in den Math. Ann. Bd. 59 dargestellt hat. Offen- 
bar kennt Verf. diese Arbeit nicht, denn sonst 
wäre ihm nicht der grundlegende Fehler unter- 
laufen, von dem noch die Rede sein wird. In 
Kap. II setzt er sich mit berechtigter Kritik mit 
den Ausführungen von Delassus, Beghin und 
Appell über die Realisierung seiner Systeme aus- 
einander. Beachtenswert erscheint der Gesichts- 
punkt, daß gefundene Integrale der Bewerungs- 
gleichungen als i. a. nichtholonome Bedingungen 
eingeführt und so von vornherein berücksichtigt 
werden können. In Kapitel III und IV setzt nun 
Verf. seine Methode auseinander und leitet die Be- 
wegungsgleichungen ab. Kapitel V handelt von 
der durch Herrn Schieldrop eingeführten deviation, 
non-holonome, das letzte Kapitel gibt noch zu- 
gehörige geometrische Betrachtungen, wobei Rie- 
mannsche Räume, die längs der Bahn in Überein- 
stimmung sind, zur Deutune des Hamiltonschen 
Prinzips eingeführt werden. (L’espace de Riemann 
de raccordement le long de la trajectoire.) 

Nun muß aber leider gesagt werden, daß die 
Methode des Verf., durch die er die älteren Unter- 
suchungen, auch eigene vereinfachen will, verfehlt 
ist. Dasselbe eilt für die Darstellung in der 
älteren Arbeit aus den Abh. der Ösloer Akademie 
1937, Nr. 11. Verf. möchte die Einführung von 
scheinbar überflüssigen Geschwindigkeitsparame- 
tern @, die durch ihr Verschwinden die nichtholo- 
nomen Bedingungen darstellen, vermeiden, und 
führt deshalb von vornherein nur h<n solche Pa- 
rameter ein, die bei Beachtung der Bedingungen 
freibleiben. Er beachtet aber nicht, daß er bei der 
Annahme döq—ödda=0, die er Seite 38 unten 
und 39 oben über die wirklichen Koordinaten trifft, 
die ö® für diejenigen ®, die selbst vermöge der 
Bedingungen verschwinden, nicht null setzen darf, 
was er aber implizite tut, da sie bei ihm gar nicht 
vorkommen. Anders ausgedrückt: Verf. möchte 
sowohl die virtuellen Verschiebungen als auch die 
Nachbarbahnen möglich haben, was bekanntlich 
ausgeschlossen ist. Die Folge ist, daß seine Glei- 
chungen 14 auf Seite 39 widerspruchsvoll sind, 
es sind n Gleichungen fürh <n Unbekannte. Wenn 
scheinbar die Auflösung dieser Gleichungen ge- 
lingt, so liegt das daran, daß er sie mit nicht be- 
stimmten Größen multipliziert, so daß das Kirgeb- 
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nis 15 unbestimmt ist. Wenn das Hauptergebnis 
auf Seite 44 wenigstens richtig gedeutet werden 
kann, indem man hinzufügt, daß in T von den 
niehtholonomen Bedingungen kein Gebrauch ge- 
macht werden darf, was man aber hineindeuten 
muß, da es aus der Ableitung nicht folgt, so liegt 
das daran, daß rein formal gesehen, gewisser- 
maßen bei Außerachtlassung der Zahl h, die Rech- 
nung in der Tat, wenn sie nämlich Sinn bekommen 
soll, so geführt wird, als wären doch n und nicht 
nur h Größen & da. Leider ist also die Ableitung 
der Gleichungen, die man dem Verf. verdankt, in 
dieser Darstellung und der genannten früheren 
unhaltbar. Wenn Verf. mit seiner Methode durch- 
kommen will, mußer auf die Übergangsgleichungen 
dögq—Ödödq=0 und damit auch auf die Lagrange- 
sche Zentralgleichung verzichten und die verall- 
gemeinerte Zentralgleichung aus den Annalen 59 
benutzen. 
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Dr. F. LOONSTRA, Analytische Unter- 
suchungen über bewertete Körper. 
101 S. Amsterdam 1941, N. V. Noord-Hollandsche 
Uitgevers Maatschappij. 

Als eine Verallgemeinerung des Rechnens mit 
absoluten Beträgen von komplexen Zahlen, wie es 
jedem Mathematiker geläufig ist, wurde durch 
Kürschak für beliebige abstrakte Körper der 
Begriff der Bewertung eingeführt (J. reine angew. 
Math., Bd. 142 |1913]). Dieser Begriff hat seitdem 
zu verschiedenen umfangreichen teils algebraischen, 
teils analytischen Untersuchungen Veranlassung 
gegeben, die auf mehrere Zeitschriften verstreut 
sind. Eine zusammenfassende Darstellung in einer 
Monographie, wie sie hier, und zwar unter Bevor- 
zugung der analytischen und Zurückstellung der 
rein algebraischen Betrachtungen gegeben wird, 
dürfte daher den Freunden und Bearbeitern dieses 
tebietes willkommen sein. 
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NACHRICHTEN 


Professor Matschoß f. 

Am 21. März dieses Jahres verstarb nach kurzer 
Krankheit im 71. Lebensjahr. Prof. Dr.-Ing. eh., Dr. 
phil. h. e. Conrad Matschoß, langjähriger 
Direktor des Vereins deutscher Ingenieure und 
Prof. für Geschichte der Technik an der Techn. 
Hochschule Berlin. Prof. Matschoß ist weiten 
Kreisen durch seine zahlreichen Arbeiten zur 
Teehnik- und Industriegeschichte bekannt gewor- 
den, von denen nur das zweibändige Werk „Ent- 
wieklung der Dampfmaschine“ (1908) und der Band 
„Männer der Technik“ (1925) genannt seien. Zur 


Förderung der Studien und zur Zusammenfassung 
der Arbeit auf diesem Gebiet gründete er 1909 die 
„Beiträge zur Geschichte der Technik und In- 


dustrie“, die seit 1989 unter dem Titel „Technik- 
geschichte“ erscheinen. 

Als Direktor des Vereins deutscher Ingenieure 
galt seine Sorge unter anderem der Nachwuchs- 
frage. Er war Mitbegründer und später Vorsitzen- 
der des Ausschusses für technisches Schulwesen 
und war ferner an der Gründung der Bewegung 
„Forschung tut not“ führend beteiligt. Der Zeit- 
schrift für angewandte Mathematik und Mechanik, 
deren erster Band 1921 im VDI-Verlag erschien, 
als Matschoß Direktor des Vereines deutscher 
Ingenieure war, brachte er stets wohlwollendes 
und verständnisvolies Interesse entgegen. 
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Prof. Dr. Jean Peters gestorben. 

Am 24. Aug. 1941 starb Prof. Dr. Jean Peters, 
Öbservator am Coppernieus -Institut in Berlin- 
Dahlem. Prof. Peters ist in weiteren Kreisen durch 
seine Tafelwerke bekannt geworden. Am 25. August 
1869 in Köln am Rhein geboren, studierte er in 
Bonn Mathematik und Astronomie. Seit 1899 ge- 
hörte er zum Astronomischen Recheninstitut in 
Berlin. Hier hat er hauptsächlich die Jahresbände 
des Berliner Astronomischen Jahrbuches bearbeitet. 
Auch der Preußische Grundkalender und das Nau- 
tische Jahrbuch sind von ihm alljährlich vorbe- 
reitet worden. Neben umfänglichen Tafeln für rein 
astronomische Aufgaben (Koordinatentransforma- 
tion, Keplerbewegung) hat Peters eine große Zahl 
Tafeln, vor allem der Kreisfunktionen, ihrer 
Logarithmen, sowie der Logarithmen der natür- 
lichen Zahlen berechnet und herausgegeben. Das 
erste dieser Tafelwerke war (zusammen mit 
J. Bauschinger) eine achtstellige Tafel 
der Logarithmen der Zahlen und der trigonometri- 
schen Funktionen (Leipzig 1910/11, Engelmann). 
Hier wird im Vorwort die Berechnung und Her- 
stellung der Tafeln eingehend beschrieben. Später 


folgten Tafeln der verschiedensten Stellenzahlen 
für die logarithmische und die maschinelle Rech- 
nung, für die Sexagesimaleinteilung, für die Stun- 
denteilung, für die Dezimalteilung des Grades, in 
der letzten Zeit auch für die Neugradteilung. Neben 
einer Tafel für das logarithmische und numerische 
Rechnen mit drei Stellen (Berlin 1913, Stankiewicz) 
steht eine Tafel der zweiundfünfzigstelligen Log- 
arithmen (Berlin 1919. Astr. Recheninstitut). Da- 
zwischen brachte er Tafeln mit 5, 6, 7, 8, 10 und 
21 Stellen heraus. Mir sind insgesamt vierzehn 
Peterssche Tafeln für das Rechnen mit Logarith- 
men und Kreisfunktionen bekannt geworden. Dazu 
kommen mehrere Multiplikationstafeln (1924, 1930), 
sowie eine Faktortafel aller Zahlen bis 100 000 (1935). 

J. Peters war ein sicherer, rascher und aus- 
dauernder Rechner. Er war groß im Organisieren 
von umfänglichen Rechenarbeiten. Seine Methoden, 
Rechen- und Druckfehler zu vermeiden, waren sehr 
gut durchdacht. Leider ist, über einen kleinen 
Kreis engster Mitarbeiter hinaus, davon kaum 
etwas bekannt geworden. Im Satzbild — auch das 
ist ein wichtiger Punkt, wenn es sich darum han- 
delt, große Zahlenmengen zu verarbeiten — sind 
die Petersschen Tafeln vorbildlich. 

Die Mehrzahl der Petersschen Tafeln ist auf 
die Bedürfnisse der rechnenden Astronomie und 
Geodäsie zugeschnitten. Gelegentlich haben aber 
auch technische Probleme den Anstoß zur Be- 
rechnung der Tafeln gegeben. Längst nicht alle 
von Peters berechneten Tafeln konnten bisher 
der Öffentlichkeit übergeben werden. So konnte 
z. B. eine achtstellige Reziprokentafel leider noch 
immer nicht gedruckt werden. 
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Preisausschreiben der Lilienthal-Gesellschaft. 

Wie in den vorhergehenden Jahren schreibt auch 
in diesem Jahre die Lilienthal-Gesellschaft für Luft- 
fahrtforschung eine Reihe von Preisaufgaben aus. 
Die Themen sind den Gebieten der Aerodynamik, 
der Flugzeugstatik, des Flugmotorenbaues (2 The- 
men), des Funkwesens und der Ausrüstung (3 The- 
men) und des Waffenwesens entnommen. Da im 
letzten Jahre trotz des Krieges eine überraschend 
große Zahl guter Arbeiten eingegangen ist, ist. 
um alle preiswürdigen Arbeiten auszeichnen zu 
können, in einzelnen Fällen die Zahl der Preise 
erhöht. Auch dieses Preisausschreiben dient dem 
großen Ziel: Verbesserung aller technischen Mittel 
der Luftwaffe. 

Themen und Bedingungen teilt das „General- 
sekretariat der Lilienthal-Gesellschaft für Luftfahrt- 
forschung, Berlin SW 11, Prinz-Albrecht-Str. 5“ mit. 
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